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ABSTRACT 



' Context. Galactic chemical evolution models are useful tools to interpret the large body of high-quality observational data on the 

' chemical composition of stars and gas in galaxies which have become available in recent years. 

, Aims. This is the second paper of a series which aims at quantifying the uncertainties in chemical evolution model predictions related 

to the underlying model assumptions. Specifically, it deals with the uncertainties due to the choice of the stellar yields. 
Methods. We adopt a widely used model for the chemical evolution of the Galaxy and test the effects of changing the stellar nucle- 
Mh' osynthesis prescriptions on the predicted evolution of several chemical species. Up-to-date results from stellar evolutionary models 

Q , are taken into account carefully. 

5_( ■ Results. We find that, except for a handful of elements whose nucleosynthesis in stars is well understood by now, large uncertainties 

still affect the model predictions. This is especially true for the majority of the iron-peak elements, but also for much more abundant 
' species such as carbon and nitrogen. The main causes of the mismatch we find among the outputs of different models assuming 

difi'erent stellar yields and among model predictions and observations are: (i) the adopted location of the mass cut in models of type II 
supernova explosions; (ii) the adopted strength and extent of hot bottom burning in models of asymptotic giant branch stars; (iii) the 
neglection of the effects of rotation on the chemical composition of the stellar surfaces; (iv) the adopted rates of mass loss and of 
^ ' (v) nuclear reactions, and (vi) the different treatments of convection. 

CO , Conclusions. Our results suggest that it is mandatory to include processes such as hot bottom burning in intermediate-mass stars 

■ and rotation in stars of all masses in accurate studies of stellar evolution and nucleosynthesis. In spite of their importance, both these 

OO ' processes still have to be better understood and characterized. As for massive stars, presupemova models computed with mass loss 

lO ' and rotation are available in the literature, but they still wait for a self-consistent coupling with the results of explosive nucleosynthesis 

I computations. 

' Key words. Galaxy: abundances - Galaxy: evolution - nuclear reactions, nucleosynthesis, abundances 

o ■ 

1. Introduction are returned to the ISM in the form of newly produced ele 

ments during the entire stella r hfetimes. F ollowing the p i oneer 



k>( . , ,. , • , , ments aurmg tne entire steiia r iireiimes. r oiiowmg tne p i oneer- 
p<; A large, everexpandmg body of high-quality observational data ■ ^^^^^ by Iben & TruranI ([1978) and Renzini & Vol? (T98ll) 
H IS available nowadays for our own Galaxy as we 1 as for a few intermediate-mass stai's (LIMSs), and Arnett ( 1 971 
. . ^ nearby galaxies, concerning the detailed chemical properties of ^ ^ai^ ([T979h for massive stars, stellar models 
their interstellar medium (ISM) and stellar popu ations. This ^^^j^ increasingly improved treatment of the relevant phys- 
wealth of data can be interpreted in the framework of galactic .^^^ processes have appeared in the literature, thus providing 
chemical evolution (GCE) models. Even though GCE is not yet a ^.^^^^^^ .^^ .^^^^ ^^^^ important differences 
well-sound asti-ophysical theory, it is still a powerful instrument ^ ^^^-^^^ ^^^^ yiel ds is the presence of mass loss 
to understand when and how a given system formed and evolved. ^^^^.^^ ^^^^j^ ^^^^^ IChiosi & Caimmil 's (fT979l) pa- 
in fact, as different chemical elements are produced on different . ^j^^ ^^^^ important work on yield s from massive stars to 
time scales in different astrophysical sites, the behaviour of key -^^^^^^ ^^^^ j^^^ ^^^^ (fT99l . Rfaitol (IT992I) 
abundance ratios for a given galactic region allows us to trace ^^^^^^ ^^^^ ^^^^ j^^^ ^^.^^^ ^ ^^^jj^^ ^.^^^ ^^^^j ^^^^^^ 
the chronology of events in that region. Abundance ratios as a ^^^^.^^ ^^^^j^ ^.^j^ ^^j^ supersolar metalUcities. In 
function of time or metallicity also provide important tests of ^j^^^ ^^^^^ ^^^.^ 1^^^ j ^^^^^-^^^ ^^^^ produced he- 
possible nucleosynthesis scenarios (Tinsley 1979). Hum and carbon. This has the effect of depressing the produc- 
^ The imnLElay^^m this game are the stellar yields, de- ^.^^ ^^^^.^^ elements, and in particular oxygen (see also 
hned by ITinsleyl (119801) as the mass fractions of stars which 
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iDrav et al.ll2003t iDrav & Toutll2003h . This strong effect is par- 
ticularly important in high metallicity environments such as the 
Galactic bulge (see McWilham et al. 2008; Cescutti et al. 2009). 
In the following years the Geneva group has produced yields 
for massive stars with a lower mass loss rate than used pre- 
viously by Maeder ( 1992), as wel l as the inclusion of stellar 
rotation (Mevnet & Maeder l2002at iHirschi et al"]l2005t iHirschil 
[2007). The most significant result was the substantial quantity 
of primary N produced by i ncluding rotation in very-low metal- 
licity massive star models ( Hirschill2005l l2007t IChiappini et al.l 
|2006|) . Unfortunately, yields including mass loss through stel- 
lar winds are available only from a limited number of stud- 
ies, whereas the majority of nucleosynthesis studies including 
a large number of chemical species do no t include mass loss 
from m assive stars. Among these we recall 'Woosley & Weaveil 
(Il995h : fUmongi & Chiefli (2003); Chiefli & Limongi (200% 
The recent study by Kobayas hi et al.l <Q.006) includes mass 
loss but no rotation. Concerning LIMSs, nucleos y nthesi s 
yields including mass loss were fro r n iRenzini & Volil (Il981h . 
Ivan den Hoek & Groenewegen (1997^, 'Marigo'('2001'). More re- 
cently Karakas & Lattanzio (2007) and Karakas (2010) have 
provided metallicity dependent yields. Concerning type la su- 
pernovae (SNela, exploding white dwarfs in binar y systems) the 



most c ommonly adopted yields are those from iNomoto et al 



(•1984", their model W7), subsequently revised bv'Iwamo to et al 
( 1999). These yields are preferred over others because they can 
very well reproduce the observed abundance pattern in SNIa 
ejecta. Empirical yields for m assive stars and SNIa were de- 
rived by iFran^ois et al.l (|2004|) . who modified the existing yields 
in order to reproduce the solar abundances and the [X/Fefl ver- 
sus [Fe/H] patterns observed in the solar vicinity. Their yields 
mainly differ from the standard ones (in this case those of 
I Woosley & Weaveiil 199 51) for Mg, Si, and the Fe-peak elements. 
However, it should be noted that no computed yields can repro- 
duce the abundance patterns of the Fe-peak elements in the solar 
vicinity. 

The yields currently available for the most abundant ele- 
ments result from updated stellar evolution and nucleosynthe- 
sis calculations, but differ quite significantly from each other 
This implies that chemical evolution models assuming differ- 
ent yields can predict quite different element abundances even 
if they assume similar galactic parameters, such as the infall and 
star formation rate (SFR) or the initial mass function (IMF). It 
is therefore important to quantify what is the effect of different 
yields on the model results. To this purpose we compare stellar 
yields currently available in the literature and describe the cor- 
responding differences in the outcome of a well tested chemical 
evolution model for our Galaxy. 

The layout of the paper is as follows. In Sect. |2] we con- 
cisely review the prescriptions about the physical ingredients of 
the stellar models adopted by different authors. This provides 
the different yield sets explored in this study. A comparison be- 
tween the different yield sets is also presented. In Sect. [3] we 
introduce the adopted GCE model. Our results are presented and 
compared with the coiTesponding observational data in Sects. |4] 
to|6l In Sect.|4]we deal with the evolution of elements from car- 
bon to zinc in the solar neighbourhood. The Galactic abundance 
gradients are briefly discussed in Sect. |5] The helium variation 



' Throughout this paper, the standard 
log (A'A/A'B)-log(A'A/yVB)o and loge(A) 



notations [A/B] = 
log{NJNn)+i2 are 
used, where Na, A'b and A'h are the number densi ties of elements A, 
B and hydrogen. Solar abundances are taken from lGrevesse & Sauvall 
unless otherwise specified. 



with metallicity is analysed in Sect. |6] Sect. [T] summarizes and 
discusses the results. 



2. Stellar yields 

In this section we provide a detailed comparison between the 
different yield sets adopted in this work. The reader is refeiTed to 
the source papers for a more thorough discussion of the adopted 
input physics. 

For LIMSs (see Sect. 12. Ill , we co nsider the yields by: (i) 
Ivan den Hoek & GroenewegenI (Il997l) . for two choices of the 
ma ss loss rate and minimum core mass for hot bottom burning; 
(ii)|Marigo'(2001), for t wo cho i ces of the mixing-length param- 
eter; (iii) Karakas & Lattanzio! (l2007h standard set, ta king the 
contrib ution from 'extra pulses' into acc ount, and (iv) KarakasI 
(2010'), who has recently recomputed Karakas & Latta nzio s 
(2007) models with updated nuclear reaction rates and enlarged 
the grid of masses and metallicities. 

Fo r massive stars (see Sect. 12.21 1. we consider the yields 
by: (i) IWooslev & Weaver! (Il995h. computed withou t mass loss 
and without rotation; (ii) Ko bayashi et alJ (l2006l) . including 
metallicity-depende nt mas s loss, and (i i i) the Geneva group 
(Mevnet & Maedeil|2002i iHirschi et al.l 120051: IHirschil 120071: 
lEkstrom et al.i2008 \ limited to the presupernova stage, but com- 

? uted with both mass loss and rotation. Older yields by iMaedeil 
1992) for solar-metallicity stars - with higher mass loss rates 
but without rotation - are considered as well. All of the adopted 
yield sets are metallicity dependent and are provided for various 
metallicities. 



2.1. Low- and intermediate-mass stars 

Stars with initial masses between ~1 and ~5-8 M© (depend- 
ing on the overshooting adopted in the stellar evolutionary mod- 
els) pass through a phase of double-shell burning, referred to 
as thermally-pulsing asymptotic giant branch (TP-AGB) phase, 
before ending their lives as white dwarfs. During the TP-AGB 
phase, they experience a number of thermal pulses followed by 
mixing episodes, the third dredge-up (TDU). In the most mas- 
sive models, hot bottom burning (HBB) can occur if the base of 
the convective envelope is exposed to temperatures hot enough 
to trigger proton-capture nucleosynthesis. During the TP-AGB 
evolution, the stars experience a very rich nucleosynthesis whose 
products (most notably He, C, N, Na, Mg and ^-process ele- 
ments) are convected to the outermost layers and eventually in- 
jected into the ISM by stellar winds. These stars are hence im- 
por tant contributors to the chemica l enrichment of galaxies. 

Ivan den Hoek & GroenewegenI (Il997l) published yields of 
*He, '^C, '^C, '^N and '^O for stars with initial masses in 
the range 0.8-8 Mq and initial compositions (Y, Z) - (0.243, 
0.001), (0.252, 0.004), (0.264, 0.008), (0.30, 0.02) and (0.34, 
0.04). They used the metall icity-depe ndent tracks provided by 
the Geneva group (e.g. iSchaller et al.lil992) to follow the stel- 
lar evolution_prior_toJheA^ with the recipes outlined 
in iGroene wegen & de Jong' ('1993') to include the effects of the 
first and second dredge-ups on the yields. Then, they followed 
the chemical evolution and mass loss up to the planetary nebula 
(PN) ejection with a synthetic TP-AGB evolution mode l. The ef- 
fect of HBB was included in an approximate way and a lReimersI 
(1975) law was assumed for the mass loss of stars on the AGB. 
These models do not consider any extra mixing. In this work 
we explore the Ivan den Hoek & GroenewegenI (1 19971) yield sets 
with mass loss scaling parameter along the AGB, ?7agb, either 



D. Romano et al.: Quantifying the uncertainties of ciiemical evolution studies. II. 



3 



vaiying with metallicity or constant. When the set with a con- 
stant ?7agb is assumed, we also investigate the effect of reducing 
the mass range in which HBB is allowed to occur 

iMarigol d2001h followed the evolution of LIMSs with ini- 
tial mass from ~0.8 to 5 M0 for three choices of the original 
composition, namely (Y, Z) = (0.240, 0.004), (0.250, 0.008) and 
(0.273, 0.019). From the zero-age main sequence up to the onset 
of the TP- AGE pha se, she has used complete stellar models by 
iGirardi et alj (l2000l) . including moderate core overshooting and 
external convection, but no extra mixing. Later phases ha ve been 
treate d with the hel p of synth etic TP-AGB models (Marig oet alJ 
11999 ). The classical FReimersI (1197 5) formula has been used to in- 
clude mass loss along the red giant branch (RGB), while mass 
loss during the TP-AGB phase has been in cluded by means of 
the sem i-empirical formalism developed bv lVassihadis & WoodI 
(Il993h . Yields of ^He, ^He, '^c, '^c, 14^^ 15^^ 16q Uq 
'^O are provided for three choices of the mixing-length param- 
eter, namely a = 1.68, 2.00 and 2.50. We investigate the effect 
of choosing the yields corresponding to either the lowest or the 

highest a value in Sect. 14.11 

The AGB yields of Kar akas & Lattanziol (l2007h and lKarakasI 

(l20Toh were calculated from detailed stellar models, where the 
structure was computed first and the nucleosynthesis calculated 
afterward using a post-processing algorithm. Yields are included 
for 77 nuclei including all stable isotopes from H to '''*S, and 
for a small group of Fe-peak nuclei. The details of this proce- 
dure and the codes used to compute the models have been pre- 
viousl y described in some detail, see for example [Karakas et al.l 
(I2009h and references therein. The models cover a range in mass 
from 1.0 Mq to 6 M© and compositions Z = 0.02, 0.008, 0.004, 
and Z = 0.0001. All models were evolved from the zero-age 
main sequence to near the tip of the TP-AGB. The TDU effi- 
ciency governs the nucleosynthesis in the lower mass models; 
this was found to vary as a function of the H-exhauste d core 
mass, metallicity, and envelope mass (see IKarakas et al.ll2002i 
for details). For example, in the Z - 0.02 models, no TDU 
was found for m < 2 Mq. For the intermediate-mass mod- 
els, the TDU was found to be efficient and the o ccuiTence of 
HBB also played a strong role in determining the final yields. 
The occurrence of HBB also depends on the initial mass and 
metallicity, with HBB occurring in lower mass stars with a de- 
crease in metallicity (at 3 Mq at Z = 10"^, whereas it only 
starts at 5 Mq at Z - 0.02). Furthermore, HBB is eventually 
shut d own by the action of mass lo ss. F or further detail s we re- 
fer to IKarakas & Lattanziol (l2007h and iKarakasI (l2010h . These 
models employ the mixing-length theory of convection with a - 
1.75. On the RG B Reimer's mass lo s s is a dopted with tjr - 
0.4. On the AGB, 'Vassiliadis & Wood' ('1993') mass loss is used 
for all stel lar mode ls in Karakas & Latta ii z i o ( 200 7) and most 
models in IKarakasI (l20ld) : in karakas! l2010r Reimer's was 
used on the AGB for models between 3 to 6 Mq with Z - 
IQ~ * (see reference for details) . The main difference between 
the [K arakas & Lattanziol (12007 ) models and those presented in 
IKarakas (.2010l) " is the choice of reaction rates used in the post- 
processing algorithm. IKarakasI (1201 Oh used an updated set of 
proton and c-capture rates that include some of the latest ex- 
perimental results for important reactions involved in the CNO 
cycl e, and the NeNa and MgAl chains. Furthermore, iKarakasI 
(I2OI O) assumed scaled -solar initial abundances for t he Z - 0.008 
and Z = 0.004 models: IKarakas & Lattanziol (l2007l) adopted ini- 
tial abundances for the Large and Small Magellanic Clouds from 
[Russell & Dopita ( 1992) which are sub-solar f or C, N, and O. 
The m ain d ifference between the yields of Kar akas & Lattanziol 
(120071) and lKarakasI ([201a) is that ~6-30 times less Na is pro- 



duced by intermedia t e-mas s models with HBB. The yields of 
IKarakas & Lattanziol (l2007l) also included the contribution from 
'extra pulses', that is thermal pulses that were not computed 
using the detailed stellar evolution code but with the help of 
a synthetic TP-AGB model. This was done because it was not 
always possible to evolve the models to the tip of the AGB 
owing to convergence difficulties. If there was sufficient enve- 
lope mass left, further t herma l pulses and TDU could poten- 
tially occur. In IKarakasI ( 1201 Ol) . the contribution from these ex- 
tra pulses was not included. This was because the synthetic 
AGB models required an assumption about the efficiency of the 
TDU for these last pulses. Full calculations have shown that the 
efficiency approaches zero as the envelope mass is decreased 
(VassiHa dis & Woodlll993t IStraniero et aT]|1997l: IKarakas et al] 
12002) . The new intermediate-mass Z = 0.0001 AGB models pre- 
sented in Karakas (2010) were evolved to small envelope masses 
(~0.1 Mq in some cases). These models did not show any de- 
crease in TDU efficiency with decreasing envelope mass, but 
in general low-me talhcity models tend to experience more effi- 
cient mixing (e.g. iBoothrovd & SackmaniJl 19881 : [Karakas et al.l 
[200l . 

2.2. Massive stars 

Only a small fraction of the stellar mass, fro m 7 to 20 per 
cent, depending on the choice of the IMF (see [Romano et al.l 
[20051 hereinafter Paper 1), is in stars with masses above 8 Mq. 
Nevertheless, these stars are major galactic polluters. They pro- 
duce almost the whole body of the a and r-process elements 
and, if rotating fast, large amounts of primary '-^ C, '"^N, ^^Ne 
and 5 -process elements at low metallicities (see iMevnet et al.l 
120091 for a review on rotating high-mass stars). Massive stars 
are also thought to be the main factories of copper in galaxies 
(Romano & Matteucci 2007). 

IWooslev & Weaveii (jl995h provided the nucleosynthetic 
yields of isotopes lighter than A = 66 (zinc) for stars of vari- 
ous metallicities (Z/Zq = 0, 0.0001, 0.01, 0.1 and 1) and masses 
in the range 1 1-40 Mq. In their models, each star is exploded 
using a piston to give a final kinetic energy of the ejecta at in- 
finity of about 10^' ergs. The explosive modifications to the pre- 
supernova nucleosynthesis are computed by taking the effects 
of neutrino irradiation into account. Mass loss is not included 
in the presupernova models. Following the onset of a success- 
ful shock, stars more massive than 30 Mq may experience con- 
siderable reimplosion of heavy elements. The amount of mat- 
ter which falls back onto the collapsed remnant depends mainly 
on the initial energy of the explosion. For stars with m > 30 
Mq, yields are given for different choices of the explosion en- 
ergy (models labelled A, B, C in .Wooslev & Weaver.1995) . We 
adopt both case A and case B, the latter corresponding to a 
slightly larger final kinetic energy of the ejecta at infinity (typ- 
ically 1.9 X 10^' ergs rather than 1.2 x lO''' ergs). We system- 
atically halve their Fe yields. This leads to a better fit to most 
of the observed abundance ratios (see Sect. [4]i. A factor of 2 of 
reduction is well within the uncertainties in modelling the ex- 
plosion. Furthermore, the reduced Fe yields are consiste nt with 
observations of several type 11 supernovae (SNell; Timm es~et al.l 
iT99l . 

It is now acknowledged that there exist two types of core- 
collapse SNe. One is normal SNell (including lb and Ic), with 
explosion energies of the order of 10^' ergs. The other is hyper- 
novae (HNe), with mo r e than 10 times larger ex plosion energies 
(IIwamotoet al.lll998h . iKobavashi et all (l2006h have calculated 
the yields for both SNell and HNe as functions of the initial 
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Table 1. Yields for massive stars from presupernova Geneva models including mass loss and rotation. Listed are: the initial metal- 
licity of the stars, Zjni, the initial stellar masses (in M0) and rotational velocities (in km s"'), mini(L'ini), and the reference papers 
from which the yields have been picked up. See the text for details. 
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Ref. 





9(500), 15(800), 25(800), 40(800), 60(800), 85(800) 


1 


10-** 


9(500), 20(600), 40(700), 60(800), 85(800) 


2 


10-5 


9(300), 15(300), 20(300), 40(300), 60(300) 


3 


0.004 


9(300), 12(300), 15(300), 20(300), 25(300), 40(300), 60(0) 


3 


0.02 


12(300), 15(300), 20(300), 25(300), 40(300), 60(300) 


4 



Referen ces. (1) Ekstrom et a U200&) ; (2) lHirschil ( l2007l) : (3) lMevnet & Maedej 
( I2002al) : (4) lHirschi et alj ( 12001 ! 



masses (from 13 to 40 Mq) and metallicities (Z - 0, 0.001, 
0.004, 0.02) of the progenitor stars. Their calculations start on 
the main sequence and proceed through core collapse, with 
the i nclusion of metallicity-dependent mass loss (de Jag er et alJ 
11988; Kudritzki et al. 1989 ). The progenitor model is exploded 
when a central density of 3 x 10'" g cm"^ is reached. For SNell, 
the final yields are obtained by setting the mass cut in such a 
way that 0.07 Mq of Fe are ejected. For HNe, the final yields are 
obtained by adjusting the free parameters involved in the 'mix- 
ing and fallback' mechanism (lUmeda & Nomoto 2002) to give 
[O/Fe] ^ 0.5 in the ejecta, as observed in extremely metal-poor 
stars. In order to include HN nucleosynthesis in GCE models, 
one should introduce a free parameter describing the fraction of 
HNe. In the present work, we analyze the two extreme cases 
where either none (shn - 0) or all (ehn = 1) stars with m > 20 
Mq explode as HNe. 

The inclusion of rotation in stellar evolutionary models con- 
siderably alters the outputs of traditional nucleosynthesis. Yields 
of ^He, '^C, '"^N and '^O are available for high-mass rotat- 
ing stars, for a fine grid of initial masses and metallicities 
(Meynet & Maeder 2002a; Hirschi et al. 2005; Hirschi 2007; 
lEkstrom et al . 2008, see Table [B- Yields of '^C, '^O and '^O 
are provided only for a subset of models. The models include 
diffusion by shears, meridional circulation and mass loss and are 
computed up to the end of the C-burning phase. The absence 
of a proper, self-consistent computation of the nucleosynthesis 
occurring at the SN stage is the major shortcoming of such mod- 
els. However, the presupernova yields of '*He and CNO nuclei 
should be almost unaffected by subse quent explosive nucleosyn- 
thesis (e.g.'Woosl ev & Weaveill995h and can thus be safely used 
in GCE models. The grid of yields for high-mass rotators that we 
use in this study has been obtained by patching together results 
presented in different papers by the same group (see Tab l efT). W e 
also test the models comp uted by [M eynet & Maedeil (l2002ah . 
iHirschi et ^ (12005 ) and E kstrom et all (120081) for non-rotating 
high-mass stars (see Sect. 14. 11 1. 



2.3. Comparison among different yield sets 

In the previous sections, we have considered a number of stel- 
lar evolution and nucleosynthesis studies. We have reviewed 
the basic input physics used in the stellar structure computa- 
tions and schematically discussed the resulting yield sets. Other 
works on stellar evolution and nucleosynthesis exist in the lit- 
erature besides the ones considered here, providing different 
grids of yields. However, these grids are computed either for 
mass/metallicity ranges smaller than the ones used here, or with- 
out the inclusion of some important physical processes affecting 
the final yields. Others are not much different from those con- 



sidered here. The latter, for ins tance, is the case of the yields 
presented by 'Izzard et al.' ('20 p4|). wh ich are very similar to the 
on es by [K arakas & Lattanzial 2007h (see figures from 2 to 7 
of iKarakas & Lattanziol 12007 ) and are thus not included in the 
present study. The Geneva group provides yields for LIMSs to 
be used to complete their grid of high-mass models. However, 
the computation of the LIMS models is stopped at the beginning 
of the TP-AGB phase. When dealing with LIMSs, the inclusion 
of the final nuclear burning phase is an essential step towards 
computing trustworthy yields of the CNO elements. Therefore, 
we prefer not to include the Geneva yields for LIMSs in our 
GCE models until they do cover the TP-AGB phase too. As a last 
example, Ventura & D' Antona (2008a b) have recently provided 
new grids of yields for intermediate-mass stars for many chem- 
ical species. With the main goal of testing the self-enrichment 
scen ario by massive AGB stars for Galactic globular clusters 
(see IVentura & D'Antonall2008al and references therein), their 
prescriptions cover a relatively narrow range of initial stellar 
masses (3-6.5 Mq) and chemical compositions (Z - 0.001 and 
0.004). As such, they are not suitable for use in chemical evolu- 
tion models of the Milky Way and other systems where stars of 
all masses and metallicity have contributed to the enrichment. 

The results from stellar nucleosynthesis are steadily improv- 
ing with time, but we are still missing a complete and homo- 
geneous set of stellar yields - the 'ideal grid' of yields com- 
puted with the same input physics for various initial masses and 
metallicities of the stars, spanning the whole range from 0.8 to 
100 Mq in initial stellar mass and from zero to (at least) twice 
solar in initial metallicity. Such an optimal yield set should also 
take into account all the relevant physical processes occurring 
in the stellar interiors and affecting the chemical composition 
of the ejected layers during the whole stellar lifetime (see iTosil 
I2OO7I for a recent appraisal of the problem). 

The potential risks of patching together different yield 
sets are highlighted in Figs. \T\ and |2] (see also Tosi 2007j). 
The left-hand panels of Fig . [T] show the solar yields of 
Ivan den Hoek & Groenewege nI (! 1997) for LIMSs combined 

1995h for massive stars (up- 



with those of Wooslev & Weaver 



per panel), the solar yields of Marigo (1200 ih for LIMSs com 
bined with those of [Maeder (1992J for He, C, N, O and 
iKobavashi et al.l (l2006h for all other e lements from m assive stars 
(middle panel) and the solar yields oflKarakasI (l2010l) for LIMSs 
combined wifli those oflHirschi et all (l2005h for He, C, N, O and 
IKobavashi et a l. (2006) for all other elements from massive stars 
(lower panel). The Marigo (2001) yield set is c ompleted with the 
homogeneous yields bv iPortinari et al.l (Il998l) for stars of initial 
mass m = 6 Mq and m - 1 Mq. Other combinations of yields 
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Fig. 1. Ejected masses of ne wly produced eleme nts as a functio n of th e stellar initial mass, from different autho rs (vdHG97 



van den Hoek & Groeneweaen 1997; WW95: Wooslev & Weaver 1995; MOl: Mari20 2001; P98: Portinari et al. 1998; M92 



Maeder.a992,; K06 :,Kobavashi et al. 2006; KL07: K arakas & Lattanzio.2007.; E08: .Ekstrom et al.. .2008; K10: lKarakal2010 ; H05 



Hirschi et al 



20051; MM02: Meynet & Maeder 2002^ and for different initial metallicities. Solid (orange) lines: He; long-dashed 



(red) lines: C; dotted (light blue) lines: N; short-dashed-long-dashed (blue) lines: O; short-dashed (magenta) lines: Na; dot-dashed 
(gray) lines: Mg. See the electronic edition of the journal for a colour version of this figure. 



for different metallicitiefl are shown in the right-hand panels of 

Fig.m 

Two problems are immediately apparent. First, a gap is 
present between the LIMS and high-mass star domains, from 
8 to 1 1 Mq or from 6 to 13 Mq, depending on the choice of the 



^ Since neither Ivan den Hoek & Groen ewegenI i 19971) nor 
iKarakas & Lattanzicl (l2007h provide yields for zero-metallicity 
stars, at Z = we adopt the yields corresponding to the low- 
est metallicities conside red by th o se aut hors, namely Z = 0.001 
for van den Hoek & Gr oeneweged ( Il997h and Z = 0.0001 for 
IKarakas & Lattanzig COOtIv 



yield sets. This gap forces us to interpolate the yields over a non- 
negligible mass interval. This mass interval is the domain of the 
so-called super- AGB stars, that have n ot reliable nucleosynthetic 
output ( Siessll2007l : iPumo et alJ 12008. and references therein). 
Moreover, as shown below, once weighted with the IMF this 
mass range is one of the major contributors to the Galactic en- 
richment. Second, several nucleosynthesis computations do not 
go beyond 40 Mq, which requires an uncertain extrapolation to 
higher masses. While, thanks to the negative slope of the IMF, 
this may be of little concern in modelling the recent chemical 
evolution of the Milky Way, it becomes extremely relevant for 
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the very early epochs, when only the most massive stars were 
polluting the ISM. 

To evaluate the contribution to the Galaxy enrichment of a 
given mass range, one must take the IMF into account. An effi- 
cient way to do that is to compute the IMF weighted yields. Let 
us first define for each element the integrated yield as: 

Pi = I mf^ (m)^(m)dm, (1) 

J mi 

where mi Mq and m„ = 100 Mq are, respectively, the lower 
and upper mass limits to the mass range of all stars formed in the 
examined region; m^'^im) is the yield (in solar masses) and (^(m) 
is the IMF. The IMF weighted yield is: 

P m"™(m)0(m)t/OT 
Wi = , (2) 



where nii and m2 are the lower and upper mass limits of the 
examined range. 

In Fig. [21 the yi elds of Fig. [1] are weighted with the IMF of 
iKroupa et al.l (Il993h . In computing the weighted yields, a linear 
interpolation is performed over the 8-11 Mq (6-13 Mq, depend- 
ing on the stellar models) mass range. When needed, the yields 
for massive stars are extrapolated to 100 Mq by keeping constant 
the last computed masses ejected as newly produced elements. 
Very often, a bump is seen in the production of specific elements 
from stars in the 8-11 (6-13) Mq interpolated mass range. In 
practice, more than 30% of all He and N newly produced by a 
stellar generation can come from stars in the 6-13 Mq interval if 
the Karakas (2010)+Geneva+Kobavashi et al. (2006) combina- 
tion of yields is adopted (Fig.|2] bottom panels). This enhanced 
contribution can be completely spurious. 
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Model 


LIMSs 


Adopted stellar yields 

Massive stars 




Comments 


1 


vdHG97, T/AGB var 


WW95, case B 




Reference model 


2 


vdHG97, T/AGB var 


WW95, case A 




Mass cut changed 


3 


vdHG97, t/agb var 


WW95, case B + M92 pre-SN yield 


s 


Winds from Z = Zq massive stars included 


4 


vdHG97, 77AGB var 


K06, fiHN = 




SNll yields changed 


5 


vdHG97, 77AGB var 


K06, £hn = 1 




HN nucleosynthesis included 


6 


vdHG97, ;7agb var 


K06, £hn = 1 + Geneva pre-SN yields, u 




Stellar rotation included 


7 


vdHG97, 77AGB var 


K06, £hn = 1 + Geneva pre-SN yields, v 


ni=0 




8 


vdHG97, 77AGB const 


K06, £hn = 1 + Geneva pre-SN yields, u 


„■ ^0 


Mass loss along the AGB changed 


9 


vdHG97, minimum HBB 


K06, £hn = 1 + Geneva pre-SN yields, u 


„■ 


HBB extent reduced 


10 


MOl, a = 1.68 


K06, £hn = 1 + Geneva pre-SN yields, u 


„i 7^0 


LIMS yields changed 


11 


MOl, a = 2.50 


K06, £hn = 1 + Geneva pre-SN yields, u 


„i ^0 


HBB strength increased 


12 


MOl, a = 1.68 


K06, £hn = 1 + [Geneva, Wini tO + M92] pre 


-SN yields 




13 


KL07, with extra pulses 


K06, £hn = 1 + Geneva pre-SN yields, u 


m 7^0 


AGB yields from detailed stellar models 


14 


KL07, with extra pulses 


K06, £hn = 1 + [Geneva, v,„, tO + M92] pre-SN yields 




15 


KIO, without extra pulses 


K06, £hn = 1 + Geneva pre-SN yields, u 


„i ^0 


Up-to-date nuclear reaction rates for LIMSs 



Note s. vdHG97:lvan den Hoek & GroenewegenI ([19971): WW9 5:IWooslev & ^yeavej ( Il995h : M92: iMaedej ( 119921) : K06:lKobavashi et alj ( 120061) ; 
MO li lMaried ( 120011) : KL07:'Karakas & L attanziol ( l2007t) : K10:'Karakas' ('2010^. The models adopting the yields bv'Marigd d200lh for LIMSs use 
self-consistent nucleosynthesis prescriptions from lPortinari et al.. ( 1998) for stars of initial mass m = 6 Mq and m-7 Mq. 



The only way to overcome these problems is getting homo- 
geneous and complete sets of yields. 

3. The chemical evolution model 

To compare the contribution to the Galaxy enrichment of the 
yields described in the previous sections, we adopt the two-infall 
model case B for the chemical evolution of the Galaxy. A thor- 
ough discussion of the a dopted formalism and basic equations 
can be found elsewhere ('Chia ppini et al.l[r997l 1200 ih . Here we 
briefly recall the overall evolutionary scenario. 

The inner halo and thick disc of the Milky Way are as- 
sumed to form on a relatively short time-scale (about 1 Gyr) 
out of a first infall episode, whe reas the thin disc forms inside- 
out ('Matteuc ci & FranpoislllQSQh on longer time-scales (7 Gyr in 
the solar vici nity and more than a Hubble time at the outermost 
radii; see also lRomano et al. "2000) during a second independent 
episode of extragalactic gas infall. The Galactic disc is approx- 
imated by several independent rings, 2 kpc wide. Radial flows 
are not considered here. Galactic outflows seem unlikely owing 
to the deep Galactic potential well, and are thus not included 
in the model. The adopted SFR is proportional to both the total 
mass and the gas surface densities. The efficiency of conversion 
of gas into stars is higher during the halo/thick-disc phase than 
during the thin-disc phase. Furthermore, it drops to zero every 
time the gas density drops below a critical density threshold. The 
instantaneous recycling approximation is relaxed, i.e. the stellar 
lifetimes are taken into account in details. As for the stellar IMF, 
the Kroupa et al. ( 1993) IMF is assumed (see Paper I, and fur- 
ther discussion in Sect. |6]l in the 0.1-100 Mo mass range. The 
rate of SNIa explosions is calculated as in Matteucci & Greggiol 
(Il986l) . SNela explode in close binary systems when a CO white 
dwarf has reached a critical mass limit because of accretion of 
hydrogen-rich matter from a main sequence or red giant com- 
panion. 

The yields for SNela are taken from llwamoto et aTl (Il999h . 
and we use their model W7. For single stars, we adopt dif- 
ferent sets of stellar yields (see Table |2] and discussion in 
Sect. |2]l. Overall, we compute 15 models, which differ only in 
the adopted nucleosynthesis prescriptions. Another couple of 
models. Models IS and lOS, have been computed with differ- 



ent assumptions on the IMF, to let the reader appreciate the 
IMF effect on the results. The y are t he same as Models 1 and 
1 0, respectively , but ass ume a IScalof s (Il986h IMF rather than 
a Kroupa et alfs (Il993h IMF. They are discussed at length in 
Sect. |6] and are not listed in Table |2] A linear interpolation of 
the yields is performed as a function of both mass and metal- 
licity. Since none of the massive star studies described in the 
previous sections includes prescriptions for a 100-Mq star (see 
Sect. 12.2b . we must extrapolate the yields available for the most 
massive stellar models. We choose to keep the masses ejected in 
form of newly produced elements fixed up to 100 Mq. 

In the following sections, we present the model predictions 
and compare them to one another and with the relevant data. 

4. Results on the abundance ratios 

In this section we discuss, element by element, the behaviour 
of several abundance ratios as a function of metallicity in the 
solar neighbourhood, [X/M] versus [M/H] . We show either iron, 
oxygen or magnesium as the metallicity indicator. 

4.1. CNO nuclei 



Carb on Carbon is known t o form through the 3-a process in 
stars (IWallerstein et al.|[l997i, and references therein). However, 
it is still debated wh ether it comes mainly from LIMSs or 



massive stars (see e.g 


. Henrv et al. 20001 Chiappini et alj 20031 


lAkerman et al.ll2004l 


Carigi et al.ll2005l; Cescutti et al.ll2009l for 



different points of view). 

In Fig. |3] we show the behaviour of [C/O] as a function of 
[O/H] (panel a) and [C/Fe] as a function of [Fe/H] (panel c) as 
traced by high-resolution observations of solar neighbourhood 
stars (symbols). The data are compared to the predictions of se- 
lected models (lines). 

In order to derive the stellar carbon abundances, spectro- 
scopists rely on high-excitation permitted C I features and on 
the forbidden [CI] line at 872.7 nm (free from non-LTE ef- 
fects), as well as on various molecular species, notably CH, 
C2 and CO. For solar-type stars, carbon abundances based 
on the forbidden [CI] line at 872.7 nm are highly accurate, 
especially if the blending weak feature affecting the line is 
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taken into account (•Bens bv & FeltzingI |2006|) . According to 
iBensbv & Fe ltzing (2006), who analyzed the forbidden [C I] line 
in a sample of 51 nearby F and G dwarf stars, [C/Fe] versus 
[Fe/H] is flat for subsolar metallicities down to [Fe/H] - -1.0, 
with no distinctive trends for thin- and thick-disc stars. At higher 
metallicities, from [Fe/H] and up to [Fe/H] ^ H-0.4, thin-disc 
members (which extend to higher metallicities than thick-disc 
members) show a shallow decline in [C/Fe]. The [C/O] versus 
[O/H] trends instead are well separated between the two discs, 
thank s to differences in the oxygen abundances. On the other 
hand. lReddv et al.l(l2006h . using permitted C I lines, find that the 
abundance ratios [C/Fe] for the thick disc are on average larger 
than the [C/Fe] for the thin disc in the metallicity range -1.0 < 
[Fe/H] < -0.4. Furthermore, [C/Fe] is found to increase with de- 
creasing [Fe/H] at sub-solar metallicities. A trend of increasing 
[C/Fe] with decreasing [Fe/H] results also from the metal-poor 
halo stars analyzed by Akerman et al. (2004). In the [C/O] versus 
[O/H] plot, an upturn in [C/O] is detecte d towards lower metal- 
licities (Fig. [3] panel a, filled triangles). lAkerman et alj (|2004|) 
corrected the oxygen abundances for non-LTE eflfects, but did 
not apply any non-LTE coiTection to the [C/O] ratios, under the 
assumption that non-LTE corrections for the permitted C I and 
OI fines are of similar size. Recent non-LTE calculations for 



43 halo stars by i Fabbian et al.l (l2009l) confirm and strengthen 
the case for an upturn in [C/O] for [O/H] < -2 (Fig. |3] panel 
a, crosses). The sparse data for extremely metal-poor halo gi- 
ants by S pite et al. (2005, Fig. |3] stars) should be given a lower 
weight, because unaccounted 3D effects l ikely aff'ect the carbon 
molecular lines employed in the analysis (lAsplundll2005h . 

In the [C/O] versus [O/H] diagram (Fig. [3] panel a), the 
predictions from Model 1 [short-d ashed (red) line], adopting 
the metallicity-dependent yields of IWooslev & WeaveJ (Il995h 
case B for massive stars, are in satisfactory agreement with the 
halo data, but fail to reproduce the solar ratio and the trend 
of increasing [C/ O] with increasing [O/H] found for the most 
metal-rich stars. iMaederf s (il992.) yields for solar-metallicity 
stars, computed with high mass loss rates but no rotation, tend 
to overestimate the [C/O] ratio at high metallicities [Model 3, 
long-dashed (bla ck) lin e]. Models 4 and 5, computed using the 
iKobavashi et aO (l2006l) yields with shn = and 1, respectively, 
underestimate the [C/O] ratio in the solar vicinity at all metallic- 
ities. Since they predict nearly the same trend, with differences 
less than 0.2 dex from one another, we show only the results 
from Model 4 [solid (blue) line]. Models including rotation in 
high-mass stars (i.e. Model 6 and Models from 8 to 15), slightly 
overestimate the [C/O] ratio in the Galactic halo, but are, in gen- 



D. Romano et al.: Quantifying the uncertainties of ciiemical evolution studies. II. 



9 



o 

^ -0.5 





-1 - 



■ir 



0.5 



O -0.5 - 



-1.5 



1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 


■ b) 


A ^ 


A 


— / / 

/ Aa 

/ / A 




a/ ^ 4^ ^ 

/a ' A 






: * ^ 

- -iz / 


/ / - 
/ J 


1 1 1 1 1 1 1/ 


/ : 
/ 1 1 



-2 -1 
[0/H] 



o 

-0.5 






Fig. 4. Carbon-to-oxygen (upper panel) and nitrogen-to-oxygen 
(lower panel) as functions of [O/H] in the solar neighbourhood. 
Data (symbols) are the same as in Fig. [3] Also shown are the pre- 
dictions from Models 1 (short-dashed lines) and 2 (solid lines). 
The two models differ only in the adopted location of the mass 
cut in models of SNII explosions. See the electronic edition of 
the journal for a colour version of this figure. 



eral, in good agreement with the thin-disc and solar data, apart 
from Models 10-12, that severely overestimate the [C/0]_ja- 
tios a t [O/H] > -0.6. These models use the yields from Marigol 
(1200 1) for LIMSs. To avoid overcrowding, in Fig. [3] we show 
the theoretical predictions only for a subset of models, namely 
Models 6 [dot-dashed (black) line], 9 [dotted (black) line] and 
10 [short-dashed-long-dashed (blue) line]. 

In the [C/Fe] versus [Fe/H] plot (Fig. [3] panel c), our mod- 
els predict a decreasing trend of [C/Fe] with increasing metal- 
licity in the metallicity domain -3.0 < [Fe/H] < -1.0. This 
trend is followed by an increase of [C/Fe] from [Fe/H] = -1.0 to 
[Fe/H] - -0.5 and, then, by either a shallow decline (Models 1, 
4, 5, 10) or a flattening of the ratio (Models 3, 6, 9, 13, 15) to- 
wards solar metallicities and up to [Fe/H] = -1-0.4. Predictions 
from Models 6, 9, 13 and 15 (the last two not shown in the fig- 
ure), in particular, a re in r emarkable agreement with data from 
iBensbv & Feltzinj (120061 filled circ l es in Fig. |3]l - free from 
non- LTE eff^ects - and lFabbian eTal] (l2009l crosses in Fig. O - 
taking appropriate non-LTE corrections into account. A notable 
exception is Model 5, computed with Kobavashi et al. (2006) 
yields for massive stars assuming that 100% of stars with m > 
20 M0 explode as HNe [shn - 1, lower (black) solid curve in 



Fig. 5. Same as Fig.H] but for Models 6 (dot-dashed lines) and 7 
(solid lines), adopting the pre-SN yields from the Geneva group 
with and without rotation, respectively. See the electronic edition 
of the journal for a colour version of this figure. 



Fig. [3] panel c], which underestimates the [C/Fe] ratio across the 
whole metallicity range. 

All in all, we find differences up to 0.5-0.6 dex in the pre- 
dicted [C/O] ratios and up to 0.8 dex in the predicted [C/Fe] 
ratios. This quantifies the uncertainties in the model predictions 
due to the uncertainties in the adopted stellar yields. The mod- 
els discussed up to now adopt the best yield sets from different 
authors, which differ quite a lot in the assumed input physics 
(see Sect. |2]i. To evaluate the effect of changing only one pa- 
rameter of stellar evolution on the predictions of GCE models, 
we run: (i) Model 2, which is the same as Model 1 but for the 
assumed final kinetic energy of the ejecta from m > 30 Mq 
stars at infinity (which is lower in Model 2 than in Model 1 ; see 
Wooslev & Weaver 1995, case A versus case B); (ii) Model 7, 
which is the same as Model 6 but for the assumed pre-SN 
yields of He and CNO elements from massive stars (which 
are now taken from the computations of the Geneva group for 
non-rotating stellar models); (iii) Model 8, which differs from 
Model 6 only in the assumed parameter of mass loss along the 
AGE (which is constant with metallicity rather than varying with 
Z in Model 8) and from Model 9 in that HBB is allowed to op- 
erate in a smaller mass range in Model 9; (iv) Model 11, which 
is the same as Model 10 but for the assumed mixing-length pa- 
rameter for intermediate-mass stars {a - 2.50 rather than 1.68 - 
this para meter mainly affects the production of ^He and '"^N due 
to HBB: iMarigoll200lh . 
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Fig. 6. Same as Fig. |4] but for Models 6 (dot-dashed lines), 8 
(solid lines) and 9 (dotted lines). Models 6 and 8 differ only in 
the adopted parameter of mass loss along the AGB, Models 8 
and 9 in the assumed mass range in which HBB can take place. 
See the electronic edition of the journal for a colour version of 
this figure. 



Fig. 7. Same as Fig.H] but for Models 10 (short-dashed-long- 
dashed lines) and 1 1 (solid lines), which differ only in the as- 
sumed mixing-length parameter for intermediate-mass star mod- 
els. See the electronic edition of the journal for a colour version 
of this figure. 



The predictions of Model 2 versus Model 1 , Model 7 versus 
Model 6, Model 8 versus Models 6 and 9, and Model 1 1 versus 
Model 10 regarding the behaviour of [C/O] versus [O/H] in the 
solar neighbourhood are shown in the top panels of Figs.|4l|5]|6] 
and|7] respectively. It is seen that changing the mass cut in mod- 
els of SNll explosions has a profound impact on the predicted 
[C/O] at all metallicities. The inclusion of rotation in models 
of high-mass stars substantially changes the predictions on the 
[C/O] ratio for [O/H] < -2. Smaller (or even negligible) ef- 
fects are seen at higher metallicity in dependence on the assumed 
mass loss and/or HBB efficiency in LlMSs. 

Nitrogen Nitrogen is synthesized th rough the CNO cycle 
in th e star hydrogen-burning layers jClavtonl 119831; lArnettI 
Il996h . It was recognized long ago that LIMSs are net produc- 
ers of (primary and secondary) nitrogen on a Galactic scale 
d^munds & Pagel 1978; Renzini & VoH 1981). While it was 
suggested years ago that metal-poor massive stars needed to sig- 
nific antly contribute to the production of primary nitrogen (see 
e.g.. lMatteucci|[T986l) . only recently has the physical mechanism 
been identified as rotation (Mevnet & Maeder 2002b). 

In the [N/Fe] versus [Fe/H] diagram, the interpretation of the 
observations and the search for any abundance trends are com- 
plicated by the fact that the stellar abundances show consider- 



able scatter at a fixed [Fe/H], which is largely attributable to the 
weakness of the usable spectral lines and/or to the T^ff sensitiv- 
ity of the derived abundances. Indeed, in late-type stars one can 
rely only on weak, high-excitation near-IR N 1 lines, also subject 
to non-LTE line formation and 3D effects. In halo stars, only 
the NH and CN molecular lines can be observed, and both have 
disadvantages. The NH band-head falls in a crowded region of 
the spectra, whereas an analysis of the CN lines requires a pri- 
ori knowledge of carbon abundances. Furthermore, both lines 
are affected by 3D effects, that for NH can result in -0.6 and 
-0.9 abundance corrections for turn-off stars at [Fe/H] = -2 and 
[Fe/H] - -3, respectivel y, with giants having not been investi- 
gated yet (Asplundl l2005L and references therein). 

The [N/O] versus [O/H] plot is another powerful diagnostic 
diagram, but once again lit erature data ne e d corr ections owing 
to non-LTE and 3D effects, ' israelian et al ] (12004 estimated ro- 
bust [N/O] ratios for 3 1 unevolved metal-poor stars using nitro- 
gen and oxygen abundances from near-UV NH and OH lines. 
Although their [N/O] (Fig. [3] panel b, empty triangles) should 
stay basica lly intact, the [O/H] values would decrease by ~ 0.5 
dex in 3D ('A splundl l2()05). More severe coiTections should af- 
fect the [ N/O] ratios deri ved by Spite et al. (2005, Fig.[3l panel 
b, stars): lAsplundl (l2005h estimates that they would decrease by 
~ 0.4 dex, while lower negative corrections of ~ 0.2 dex should 
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apply to the [O/H] values. This would lead to a slowly increasing 
trend of [N/O] with [O/H] even at the lowest metallicities. 

Being aware of all the shortcomings affecting both the data 
and the models, we show the predictions of some of our models 
for [N/O] versus [O/H] and [N/Fe] versus [Fe/H] in Fig.[3j pan- 
els b and d, respectively. Other predictions from models where 
the parameters of stellar evolution are changed one at a time, are 
shown in Figs.|4]to|7](only for [N/ O] versus [O/H]). Model pre- 
diction s are compared to data from Spite et aP (l2005h. iLai et al.l 
(l2008h . Ilsraelian et al.l (l2004 . and lReddv etal] (l2003h . It is ap- 
parent that only Models 6 to 15, adopting the N yields from 
low-metallicity massive stars from the Geneva group, are able to 
explain the current measuremen ts of [N/O] and [N/F e] in very 
metal-poor halo stars (see also lEkstrom et al.l l2008l their fig- 
ure 9). The most massive zero-metallicity stars produce great 
amounts o f N iiidependently of the assumed stellar rotation (see 
lEkstrom et al.ll20()8i their table 4). As a consequence, a plateau 
in [N/O], [N/O] ~ -0.9, is found for [O/H] < -2.5 when as- 
suming the Geneva group yields (Fig.|5] lower panel), indepen- 
dently of whether zero-metallicity stars rotate fast or do not ro- 
tate at all. However, all the models severely underestimate the 
[N/O] ([N/Fe]) ratio in the metalHcity range -2 < [O/H] < -1 
(-3 < [Fe/H] < -1). A contribution to the N enrichment from 
the winds of super- AGB stars could improve the model predic- 
tions. On the other hand, in the framework of our models, a 
range in [Fe/H] from -3 to -1 coiTesponds to a range in overall 
metallicity Z from about 4x10"^ to about 0.002. In this metal- 
licity range, we are using the yields obtained by interpolating 
the grids computed by Me vnet & Ma eder (2002a) for massive 
slow rotators with Z - 10"^ and Z - 0.004. Thus, the predicted 
deficienc y of nitrogen could u nravel th e need for an upward re- 
vision of iMevnet & Maedeit s (l2002ah N yields. This could be 
achieved by making the stars rotate faster at those metallicities 
(the faster the stellar rot ation, the larger the primary N produc- 
tion ; |MiineL&Maide3^2b ) . 

Changing the prescriptions about N production in LIMSs has 
a profound impact on the theoretical predictions at higher metal- 
Ucities, [O/H] > -1.4 and [Fe/H] > -2.0. It is clear that both 
changing the rate of mass loss along the AGB and the efliciency 
of HBB in models of intermediate-mass stars play a significant 
role in shaping the model predictions for N (see the lower pan- 
els of Figs. |6] and |7] respectively). Differences up to 0.7 dex in 
both [N/O] and [N/Fe] are found at intermediate metallicities for 
models keeping the same nucleosynthesis prescriptions for mas- 
sive stars and changing those for LIMSs (cfr. predictions from 
Models 9 and 10, dotted and short-dashed-long-dashed lines, re- 
spectively, in Fig. [3] par iels b and d). M odel 15, which assumes 
the revised N yields bv lKaraka"sl (1201 0) for LIMSs, is the only 
one which fits best both the solar [N/O] and [N/Fe] ratios (Fig.|3] 
thick solid line). 



Oxygen Next to hydrogen and helium, oxygen is the third most 
abundant element in the universe. It forms essentially in the in- 
teriors of ma ssive stars, through hydrostatic burning of He, C, 
and Ne (e.g. lArne till 1996h . Available indicators of oxygen abun- 
dances in stars include the OI triplet at ~7774 A (mainly ob- 
served in FG dwarfs and subgiants), the forbidden [OI] lines 
at 6300 and 6363 A (detected in giants and cool subgiants, and 
barely detectable in dwarfs), and the OH lines in the UV (ob- 
served in FGK metal-poor stars) and in the IR (only observed 
in cool stars with Teff < 5000 K). The features at 7774 A are 
strong and located in a clean part of the spectrum but subject 
to large and poorly understood deviations from LTE conditions. 
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Fig. 8. [O/Fe] versus [Fe/H] in the solar neighbourhood. 
Theoretical predictions are from Models 1 [short-dashed (red) 
lines], 3 [long-dashed (black) lines], 4 [upper (blue) solid 
lines], 5 [lower (black) solid lines] and 6 [dot-dashed (black) 
lines]. Upper p anel: oxygen abun dances from the forbidden 



[OI] lines from ICavrel et all (12004 stars). lOratton et al.l (I2003L 
upside-down triangl es), iReddv et al.l (12003 , open squares) and 



iBensbv et al.l ('2005, filled circles). Lower panel: oxygen abun- 
dances f rom other ind ic ators from lLai et al. (2008, open pen- 
tagons), iFabbian et aH ( l2009l crosse s), [Gratton et al. (200|, 
ups ide-down trian gles), Reddv et^ (|2006|, fill ed squares), 
Red dv et alJ (I2003L open squa res), iRamfrez et aP ([2007, trian- 
gles) and iBensbv et al.l (2005, filled circles). See the electronic 
edition of the journal for a colour version of this figure. 



The forbidden lines are free from non-LTE effects but weak 
and blended by lines from other chemical species. Furthermore, 
switching from ID to 3D models changes the [O I] -based abun- 
dances by as much as -0.2 dex in turn-off stars with [Fe/H] = 
-2; a similar behaviour can be anticipated for metal-poor giants 
dNissen et al. 112002.). The OH lines are difficult to observe and, 
in the UV domain, are also affected by surface inhomogeneities. 
The most reliable indicators of oxygen abundances in stars seem 
to be the forbidden [O I] lines (see volume 45 issue 8 of New 
Astronomy), especially if spectra of high resolution and high 
S/N are acquired a nd analyzed using accurate atomic line data 
dNissen et al ] l200l . 

In Fig. [8] we compare some of our oxygen model predictions 
to observations of the forbidden [O I] line as well as of other oxy- 
gen lines (panels a and b, respectively), concentrating on models 
which adopt the best yield sets suggested by different authors. 
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Fig. 9. Upper panel: predictions on [O/Fe] versus [Fe/H] from 
Models 1 and 2 [short-dashed (red) and solid (black) lines, re- 
spectively] are compared one to another and with data from the 
forbidden [OI] lines (see Fig. [8] caption for references). Lower 
panel: same as the upper panel, but for Models 6 and 7 [dot- 
dashed and solid (black) lines, respectively]. See the electronic 
edition of the journal for a colour version of this figure. 

Models 1, 6, 8, 9, 10, 11, 13 and 150 forecast an [O/Fe] ratio 
which, starting from [O/Fe] ^ -i-l at [Fe/H] ^ -4, slopes gently 
down to subsolar values at the highest metallicities. This is in ex- 
cellent agreement with the observations, especially if the [OI]- 
based abundances are considered for comparison. The gradual 
increase of [O/Fe] with decreasing [Fe/H] in the [Fe/H] < -1 
regime stems from the larger oxygen-to-iron ratios in the ejecta 
of more massive, more metal-poor SNell. The quick decrease 
of [O/Fe] above [Fe/H] ^ -1 is due instead to the overwhelm- 
ing contribution to Jhe Fe production by SNela. Models 3, 12 
and 14, assuming IMaedeif s (119921) yields for solar-metallicity 
massive stars, predict a decrease of the [O/Fe] ratio at high 
meta lhcities definitely s teeper than indicated by trustworthy data 
from Bensbv et al.] (12005, Fig. [8] panel a, filled circles). In or- 
der to avoid overcrowding, only the predictions from Model 3 
are shown [Fig.Hl long-dashed (black) cu rves]. Model 4, adopt- 
ing the yields by iKobavashi et al.l (|2006|) for normal SNe, fits 
the disc and Sun data nicely [Fig. [H upper (blue) solid curve]. 
However, it slightly overestimates the halo data. Model 5, in- 
stead, assuming that all stars above 20 M© explode as HNe with 



^ For sake of clarity, only the predictions from Models 1 and 6 are 
shown in Fig. [8] [short-dashed (red) and dot-dashed (black) lines, re- 
spectively]. 



much higher explosion energies, predicts an almost flat [O/Fe] 
versus [Fe/H] trend in the early halo [Fig. |8] lower (black) solid 
line]. Though most of the observed [O/Fe] ratios lie above the 
theoretical predictions, there are some stars whose abundances 
are consistent with a scenario of formation from pure HN ejecta. 

It is common belief that the synthesis of O in stars is well 
understood. Yet, by using different prescriptions on the O pro- 
duction in stars, we find that the predicted [O/Fe] ratio may vary 
by as much as 0.6 dex at the lowest metallicities. The difference 
reduces to 0.2 dex at solar metallicity. At the origin of the un- 
certainties in the predictions of our GCE model there is, once 
again, the complex interplay among the different assumptions 
made by different authors working on stellar evolution and nu- 
cleosynthesis to treat the different processes regulating the evo- 
lution of the stars. In Fig.|9] upper panel, we compare the predic- 
tions of Models 1 and 2, which differ only in the location of the 
mass cut in models of SNll explosions. It is seen that a variation 
in this parameter may have a huge impact (up to 0.5 dex) on the 
predicted [O/Fe] at the lowest metallicities. Stellar rotation has 
a lower - still non-negligible - impact (see Fig.|9] lower panel). 



4.2. Odd-Z elements 

Sodium Sodium is an odd-Z element, synthesized mostly dur- 
ing hydrostatic carbon burning, and partly during hydrogen 
burning throu gh the NeNa cycle; the .s-process produces some 
Na as well (ClavtonI l2003h . Adopting the LTE approximation 
in standard abundance analyses of Na can lead to substantial 
overestimates of the Na abundances. Non-LTE abundance cor- 
rections can be very large, typically -0. 1 dex at solar metallic- 
ity and -0.5 dex or more at [Fe/H] around -2, for the 589 and 
819 nm doublets; the departures from LTE become less severe 
for [ Fe/H] < -2 as the lin e formation shifts to deeper layers 
(e.g. iBaumiiller et"al]|1998l) . The weaker 568 and 615 nm dou- 
blets are less sensitive to non-LTE effects. Departures from LTE 
are most likely the culprits for the discordant values of [ Na/Fe] 
reported in the literature for metal-poor stars (see .Reddv et al.l 
2006, and references therein). From their LTE analysis of the 
Nal lines at 6154 and 6160 A, 'Reddy et al."(2006) find a hint for 
an increase of [Na/Fe] with increasing [Fe/H] in the metallicity 
range -1.0 < [Fe/H] < -0.6, followed by a decrease towards so- 
lar values for -0.6 < [Fe/H] < 0. Such a trend is confirmed by 
Shi et al. (2004) usin g non-LTE sta t istical equilibrium calcula- 
tions. For [Fe/H] > 0. iBensbv et"aD (l2003l '2005) observ e a rise 
in the [Na/Fe] trend, which was also seen by Edvardsson et al.l 
(Il993h . As for the lowest metallicities, when only extremely 
metal-poor dwarfs and unmixed giants are considered and the 
non-LTE corrections of the Na abundances are taken into ac- 
count, one gets a flat trend of the [Na/Fe] ratio in the metallicity 
regime -4.0 < [Fe/H] < -2.5, [Na/Fe]NLTE = -0.21 +0.13 dex, 
and a small sta r-to-star variation, comp arable to that found for 
other elements (lAndrievskv et al.ll2007i) . 

As one can see in Fig.[10l panel a, only Mo del 4 [upper solid 
(blue) curve], corresponding to the yields by IKobavashi et al.l 
(2006) for massive stars without HNe, predicts a [Na/Fe] ratio 
in agreement with the halo data. Models including nucleosyn- 
thesis from HNe, namely. Models from 5 to 15 [for sake of clar- 
ity, only the predictions from Model 5 are shown in the figure, 
lower solid (black) curve], underestimate the [Na/Fe] ratio in the 
halo. Mod el 1 [short-dashed (red) line], adopting the yields of 
IWooslev & Weaver ( 1995) case B, provides an acceptable fit to 
the halo ratios, but only in the -2.5 < [Fe/H] < -1 metallicity 
range. All in all, at the lowest metallicities we find differences up 
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to 1 .4 dex in the predicted [Na/Fe] ratio which, similarly to what 
we found for the [O/Fe] ratio, are largely driven by the adopted 
mass cut and energy of the explosion in models of massive stars. 
These differences reduce to ~0.2 dex at [Fe/H] > -2. 

Models from 1 to 12 do not account for Na synthesis from 
LIMSs. If Na production from LIMSs is considered, very differ- 
ent predictions for Na evolution in the disc are found, d ependin 
on the adopted yields. Model 15, adopting the revised ^araka; 
(120101) yields of Na from LIMSs (Fig. [ID panel a, thick solid 
line) predicts an almost negligible production of sodium from 
the se stellar sources on a Gala ctic scale. The adoption of the 
old iKarakas & Lattanzid (l2007h yields, instead, leads to a huge 
Na production from LIMSs in the disc (Fig. (TT] panel a, thick 
dotted line). These differences in model predictions for Na are 
es sentially due to the up-to-date nuclear reaction rates adopted 
bv lKarakasI (1201 Ol see discussion in Sect. |2l). None of the mod- 
els, however, is able to explain the intriguing upward trend of 
[Na/Fe ] at [Fe/H] > 0, first pointe d out by [Edvardss on et alJ 
(Il993h and clearly seen in the Bensby_eLaL| (f2005) data (filled 
circles in Fig. [TTT l. None of the models includes Na production 
from rotating massive stars. We will come back to this point in 
Sect. 12] 



Aluminium Aluminiu m is a product of hydro static carbon 
and neon burning (e.g. IWooslev & Weaver! [l995h . It can also 
be produced during hydrogen burning in the MgAl chain. 
Observationally, Al abundances need large non-LTE corrections: 
at solar [Fe/H], they amount to approximately +QA dex for the 
Al I 396 nm resonance line, but reach -1-0.4 dex or even more 
for [Fe/H] < -1.0; the maximum correction, -1-0.8 dex, is re- 
quired for mildly metal-poor turn-off stars in which the line is 



varied the rate of the nuclear reaction ^'*Mg(p, 7)^^A1 by orders 
of magnitude outside of the published uncertainties. 



saturated (Asplund 20 05|)- For [Fe/H] > -1 .0, Al shows a 'Mg 
like behaviour' ( Bens bv et all I2003L l200l iReddv et ; 



iws a_Mg- 
"an i2006h . 



At very low metallicities, [Fe/H] < -2.8, the relation [Al/Fe] 
versus [Fe/H] bec omes rather flat, [Al/Fe]NLTE = -0.06 + 0.10 
(lAndrievsk v et al.' 2008): the scatter around the mean value is 
small er than the scatter of the LTE determinations dCavrel et al.l 
I2OO4 I) and can be explained by measurement errors only. A few 
sparse data are available in the intermediate metallicity regime. 
Here, the non-LTE [Al/Mg] ratio displays an intriguing step- 
like difference between stars belonging to the Galactic halo and 
disc populations, that might prove useful as a fu rther member- 
ship diagnostic criterion (see lGehren et al.l[2006l and discussion 
therein). We will come back to this issue later. 

In the low -metallicity dom a in, the best-fitting models are 
those adopting 'Kobav ashi et al.l ('2006') yields for massive stars 
with £hn - 1 (see Table |2j only predictions from Model 5 
are shown in Fig. [TOl i. Such models predict an almost constant 
[Al/Fe] ratio during the first 370 Myr of Galactic evolution. 
Then, the ratio increases mildly, owing to the metal dependency 
of Al production in massive stars, and flattens again around a 
Galactic age of 1 Gyr, because of the increasing Fe production 
from SNela. None of the models can explain the behaviour of 
[Al/Fe] in the Galactic disc. The activation of the MgAl chain 
in the central regions of rotating massive stars on the main se- 
quence, followed by transportation of Al-rich matter to the outer 
envelope and ejection of the outermost layers by slow stellar 
winds, might lead to a substa ntial modification of the Al yields 
used here (D ecressin et al.l2007i) . However, complete grids of Al 
yields from rotating massive stars are not available yet. Also, one 
must be aware that, in order to reproduce the Mg-Al anticorre- 
lation observed in globular clusters stars. iDecressin et al.l(l2007h 
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Fig. 10. [X/Fe] versus [Fe/H] relations for sodium (panel a), 
aluminium (panel b) and potassium (panel c) in the solar neigh- 
bourhood. Theoretical predictions are from Models 1 [short- 
dashed (red) lines], 4 [upper solid ( blue) lines] and 5 [lower 
solid (black) lines]. Data are from lAndrievskv et al. | ( 20071 
'20081 ^ig stars, for Na, Al and K respe ctively). iLai et all 
(20081 open pentagons). iGratton et"al1 (|2003[ upside-down tri- 
angles), Gehren et al. (2006, open triangles), 'zhang et al. (20061 
open circles), Reddy et al, (2006, filled square s), Rcddv et aP 
(2003t open squares) and [B ensbv et il] ( I2005L filled circles). 
LTE abundances from 'Lai et al. (2O08Hiave been corrected by 
-0.5 dex for Na (fe aumuller et al.l ll998l) and -1-0.65 dex for Al 
dBaumiiller & Gehrenl 1 1997ft . See the electronic edition of the 
journal for a colour version of this figure. 
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Fig. 11. Same as Fig. [TO] panels a and b, but for Models 5 (thin 
solid lines), 13 (thick dotted lines) and 15 (thick solid lines). 
The predictions for [Al/Fe] versus [Fe/H] from Model 15 lay on 
those from Model 13. See the electronic edition of the journal 
for a colour version of this figure. 



The role of Al production from LIMSs should be b etter 
investigated as well: the recent study by Karakasi (1201 Oh as- 
cribes almost no Al production to LIMSs (Fig.[TTl lower panel, 
thick solid line), confirmin g results from previous work by 
iKarakas & Lattanziol (l2007l Fig. [TT] lower panel, thick dotted 
line). If AGB stars are responsible for the globular cluster Mg-Al 
anticorrelation, then it is likely that AGB stars as a whole synthe- 
size more Al than predicted by current models dMarcolini et alJ 
|2009>) . Should a significant amount of Al come from LIMSs, a 
sudden rise of the [Al/Fe] ratio in field Galactic stars would be 
obtained starting at [Fe/H] ~ -1.5 (cfr. Fig.[TT] panel a, thick 
dotted line, showing the results of Model 13, with a great Na 
production from LIMSs), in agreement with the trend inferred 
from the observations. 

Although the uncertainties in the GCE model predictions for 
Al are by far less dramatic than those for Na, especially in the 
high metallicity domain, the comparison with the available Al 
abundance data clearly shows that the synthesis of Al at disc 
metallicities is still far from having been understood. 

Potassium The alkali element potassium is mostly produced by 
a combination of hydrostatic oxygen shell burning and explosive 
oxygen burnin g in proportions that var y depending on the stel- 
lar mass (e.g. lWooslev & Weaverll 19951 and references therein). 
LTE abundances obtained from the KI 769.9 nm resonance line 



Fig. 12. [Mg/Fe] versus [Fe/H] (top panel) and [Fe/Mg] versus 
[Mg/H] (bottom panel) in the solar neighbourhood. Theoretical 
predictions are from Models 1 [short-dashed (red) lines], 4 [up- 
per solid (blue) line in panel a and lower solid (blue) line in 
panel b] and 5 [lower solid (black) line in p anel a and upper 
solid (black) line in pa n el b]. Data are from [Andrievskv et al 
(2010, stars), iLai et alj (l2008l open pentago ns), Gratton et al 
(120031 upside-down triang les), iGehren et al. (2006, op en tri- 
angles), iReddv et al.l (l2006l filled squar es), Reddy et alj (I2003L 
open squares) and iBensbv et al.l ( I2005L filled circles). See the 
electronic edition of the journal for a colour version of this fig- 
ure. 



are overestimated by 0.4-0.7 dex over a wide range of metal- 
licities; non-LTE corrections are smaller, but remain significant 
(from -0.15 dex to -0.3 dex), for the AA 12522, 12432, an d 
1 1769 A weak subordinate lines ( llvanova & Shiman skii 2000). 

In Fig.[TOl panel c, we compare our potassium model predic- 
tion s to the dat a by And rievskv et all ( l2010h . lZhang et al.l (l2006h 
andl Reddv et al.l (12003). Non-LTE corrections are taken into ac- 
count in the first two studies. Variations up to 1 dex in [K/Fe] 
are found among different model predictions. The comparison 
between data and model predictions shows the absolute inade- 
quacy of the adopted stellar yields. 

4.3. Alpha elements 

Magnesium Most of magnesium comes from hydrostatic car- 
b on burning and explosiv e neon burning in massive stars (e.g. 
IWooslev & Weaver! 1 19951) . Stellar spectroscopists may rely on 
several Mgl lines, as well as on a few Mgll lines, to derive 



D. Romano et al.: Quantifying the uncertainties of ciiemical evolution studies. II. 



15 



magnesium abundances. Non-LTE corrections are of little con- 
cern in this case; they are negative for the high-excitation Mg II 
lines and the Mgl b triplet at 517 nm (-0.05 dex and -0.2 dex 
at most at low [Fe/H], respectively), while they are positive for 
the remaining lines, the maximum correction being +0.2 dex 
at low rFe/H1 (Asplun3|200l and references therein; but see 
lAndrievskv et al.l2010L pointing to somehow larger corrections). 

GCE models adop ting the yields of Mg by 
IWooslev & Weaveii (Il995h for massive stars underestimate 
the abundance of magnesium in the discs of the Galaxy (see 
Fig. [12] sho rt-dashed lines, referring to ca se B; see also 
fTimmes et afl[T995l lG oswami & Prantzos 2000), unless ad hoc 
adjustments are made to the yields (Fr ancois et aL.,2004.) . While 
detailed nucleosynthesis calculations s till have to e s tablish 
whether the modifications proposed by i Franpois et al.l (|2004|) 
have a physical ground, HNe can provide an alternate solution: 




Fig. 13. [X/Mg] versus [Fe/H] relations for sodium (panel 
a) and aluminium (panel b) in the solar neighbourhood. 
Theoretical predictions are from Models 1 [short-dashed (red) 
lines], 4 [upper solid (blue) lines], 5 [lower solid (black) 
lines], 13 (th ick dotted line) and 15 (thick so lid lines). 
Data are from lAndrievskv etatl (l2007l l2008i l20ia st a rs, for 
Na, Al and Mg abundances, respec tively), iLai et a D (l2008l 
open pentago ns), 'Grat ton et al.l ('2003', u pside-dow n triangles), 
iGehren et aP (2006, open triangles), R eddv et a l.l (^2006, filled 
square s), Reddv et al. (2003, open squares) and Bensbv et al. 
(I2OO5L filled circles). LTE abundances from Lai et al.i (|2008j) 
have been corrected b y -0.5 dex for N a (Baumiille r et alJI 19981) 
and +0.65 dex for Al ( Baumiiller & Gehrenll 19971) . See the elec- 
tronic edition of the journal for a colour version of this figure. 



iKobavashi et alJ (l2006h have obtained a nice fit to the [Mg/Fe] 
versus [Fe/H] relation traced by solar neighbourhood stars 
using a GCE model for the solar vicinity which adopts their 
new nucleosynthetic yields for SNell and HNe, shn = 0.5 
for the HN fraction and a Salpeter IMF. Yet, their model can 
not reproduce val ues of [Mg/H] exceedi ng solar, at odds with 
observations (see Kobavas hi et al.l 12006'. their figure 9). Here 
we reasses the problem of Mg evo lution, adopting a d ifferent 
GCE model, with a steeper IMF dKroupa et alj Il993l) . better 
suited to describe the stellar populations of the Milky Way field 
(see Paper I, and references therein), and slightly different mass 
limits for the IMF, m/ = 0.1 M0 and ni u - 100 Mp), rather tha n 
mi = 0.07 M0 and m„ = 50 M© as in IKobavashi et a ll (l2006l). 
If we trust the non-LTE Mg abundances recently computed by 
lAndrievskv et al.l ( 1201 Ol) for the most metal-poor stars, the best 
fit to the observed [Mg/Fe] vers us [Fe/HI and [F e/Mg] versus 
[Mg/H] relations is obtained with'Ko bavashi et al] (|2006) yields 
and Shn values intermediate between the two extreme cases 
analysed here (solid lines in Fig.[T2ll. On the other hand, if the 
LTE Mg abundances by 'Lai et alj (l2008l) are the ones to be 
trusted, the extreme case of shn = 1 is favoured instead. None 
of the models, however, is able to reproduce the highest [Mg/H] 
values and the flattening of the relations at nearly solar/higher 
than solar metallicities. This could suggest the need for either 
a revision of current SNII and/or HN yields for solar and/or 
higher than solar metallicity stars, or larger contributions to Mg 
production from SNela, or significant Mg synthesis in LIMSs, 
or a combination of all these factors. Clearly, it is not safe to use 
Mg as the prime metallicity indicator in metal-rich galaxies. 

It has been recently pointed out by Ge hren et alj (2006) that, 
while the trend of [Na/Mg] versus [Fe/H] is witnessing a history 
of continuous enrichment from the halo to the discs, the [Al/Mg] 
ratio displays a much more interesting step-like behaviour, with 

halo stars sharing a common value of ([Al/Mg]) 0.45 and the 

two disc populations mostly consistent with ([Al/Mg]) — hO.lO. 
At the lowest metallicities, say -3.6 < [Fe/H] < -2.5, ex- 

tremely metal-poor stars share a lower value, ([Al/Mg]) 0.7 

(lAndrievskv etal.1 l2010i) . In Fig. \13\ we show these data, to- 
gether with measurements from other authors. Superimposed are 
our theoretical tracks, for a significant subset of models (lines). 
Three results are immediately apparent from this figure. The 
first is that all models overpredict the [Na/Mg] and [Al/Mg] ra- 
tios at [Fe/H] > 0, thus reinforcing the demand for a Mg pro- 
duction higher than predicted by current GCE models at late 
times. The second important result is that t he step-like behaviou r 
of [Al/Mg] versus [Fe/He], suggested bv iGehren et all (l2006l) . 
can not be reproduced by the models. A jump in [Al/Mg] at 
[Fe/H] X -1 could be achieved if LIMSs were efficient Al pro- 
ducers (cfr. Model 13 predictions for Na - thick dotted line in 
Fig. [13] panel a). Finally, it is seen that different models predict 
almost the same [Al/Mg] versus [Fe/H] behaviour in the solar 
neighbourhood (Fig. [131 bottom panel), whereas the theoretical 
[Na/Mg] versus [Fe/H] relations may differ by as much as 0.5- 
1.0 dex or even more at the lowest metallicities (Fig. [13] top 
panel). The best description of the observational [Na/Mg] ver- 
sus [Fe/H] diagram would be obtained with a model assuming 
a set of yields for massive stars int ermediate between the two 
sets used in Models 1 and 5, namely , IWooslev & Weaveil (Il995h 
case B and IKobavashi et all (I2OO6I) with Shn - 1, respectively 
(see TableO. 



Silicon Silicon is produced mostly during oxygen burning and 
its final abundance in the stellar ejecta is sensitive to a variety of 
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Fig. 14. Predicted behaviour of the relative ratios of several a elements to iron (panel a: silicon, panel b: sulphur, panel c: calcium, 
panel d: titanium) with respect to the relative iron abundance in the solar neighbor hood, for Mo dels 1 [short -dashed (r ed) lines], 4 

(12008. open 



[upper solid (blue) lines] and 5 [ lower soli d (black) lines ] . Also shown are data from lCavrel et al.l(,2004. star s). iLai et al 



pentag ons), iNissen et all (120071 cro sses). iGratton et~al] (l2003l upside-down triang les), ICaffau et al.l (l2005l triangles). iReddv et al.l 
(2006, filled squares'). iReddv et all (120031 open squares') and Bensbv et al.l (l2005l filled circles). See the electronic edition of the 
journal for a colour version of this figure. 



factor s regulating the evolution and explosion of the progenitor 
stars (I'Wooslev & "WeaverlllQQSh . Observationally, silicon abun- 
dances are aifected by contamination from CH and H6 lines; de- 
partures from LT E do not seem t o be a major concern, at least 
for Sun-like stars dShi et al.ll2008l) . 

At the lowest metallicities, where the signature of low-Z 
SNell is clearly seen, the [Si/Fe] ratio predicted by our GCE 
mod el may vary by up to 1 dex . Models 1 and 5, compute d 
with Wooslev & 'Weaveit s (Il995h case B and lKobavashi et all 's 
shn = 1 yields for massive stars, respectively, predict a 
mildly declining trend of [Si/Fe] versus [Fe/H] over the whole 
metallicity interval covered by the observations [Fig. [14] panel 
a, short-dashed (red) and lower solid (black) lines, respectively]. 
The theoretical trend is in good agreement with the available 
data, apart from th e flattening traced b y observations at super- 
solar metallicities (iBensbv et al.l 12005 *). which might require a 
revision of current stellar yields from high-metallicity massive 
stars and/or SNela. A very steep decline of the [Si/Fe] ratio with 
metallicity in the -4 < [Fe/H] < - 1 metallicity rang e is obt ained 
instead by assuming the yields of iKobavashi et all (l2006h with 
£hn = [Model 4, Fig. [14] panel a, upper solid (blue) line], at 
odds with observations. This would indicate that in the very early 
halo phases most massive stars must explode as HNe rather than 
normal SNell. 



Sulphur Like silicon, sulphur is an intermediate-mass element 
whose production is directly attributable to oxygen burning, ei- 
ther at the center or in the convective shel ls of massive stars, 
or during explosive nucleosynthesis (e.g. W ooslev & Weaver! 
[19951 iLimongj&Chieffi 2003). As a volatile element, it is 
hardly aff'ected by depletion onto dust particles, which makes 
it extremely useful in deciphering the ch emical enrich ment of 
damped Lyman a systems (see e.g. Nisser retal1l2004l) . Hence, 
unraveling the history of S enrichment in the early Galaxy has 
become a fundamental step towards our understanding of the 
high-redshift universe. However, sulphur abundan ce determina- 
tions in Galactic stars are currently highly debated. 'Takeda et alj 
(2005) find non-LTE abundance corrections to be negligible for 
the AA 8693, 8694 A S 1 Unes of multiplet 6, while the AA 9212, 
9228, 9237 A S 1 lines of multiplet 1 suffer significant nega- 
tive non-LTE corrections (up to 0.2-0.3 dex). In the very metal- 
poo r metalhcity re gime, a marked discordance is reported by 
Takeda et al.l (12005 ) between the [S/Fe] values from the two 
abundance indicators: while the lines of multiplet 1 attain a 
nearly flat plateau (or even a slight downward bending), the lines 
of multiplet 6 show an ever-increasing trend with decreasing 
metallicity. A more recent analysis by Niss en et al.l (120071) . how- 
ever, based on high S/N ratio spectra of 40 metal-poor halo stars, 
using both weak (A 8694.6) and strong (AA 9212.9, 9237.5) S I 
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lines, and taking non-LTE corrections into account, finds a flat 
trend of [S/Fe] versus [Fe/H] and no stars with [S/Fe] in excess 
of 0.60 dex. 

In Fig. [ 1 4] pan el b, the sulphur abundances measured by 
iNissen et aL r?2007, crosses ) are com pared to the much more 
dispersed measurements by ICaff^au et al. (2005, triangles) from 
multiplets 1, 6 and 8. Disc data at higher metallicities are from 
iReddv et al.l (|2003 [ open squares). Also shown in Fig. 1141 panel 
b, are the predictions of our GCE models for the evolution of 
[S/Fe] versus [Fe/H] (s olid and dashed line s ). M odels 1 and 5, 
adopting the yie l ds of IWooslev & Weaver! (fT995.) case B and 
iKobavashi etal] (l2006h with e = 1 for massive stars, respec- 
tively, favour a nearly flat trend of [S/Fe] versus [Fe/H] in the 
metallicity regim e -3 < [Fe/H] < - 1 . M odel 4, instead, adopt- 
ing the yields of Kobavashi et al.l (|2006|) with s — for mas- 
sive stars, predicts a steeply decreasing trend of [S/Fe] versus 
[Fe/H] in the same metallicity range. It is worth noticing that 
the outcome from Model 4 is at variance with observations if 
the data by Nissen et al. (2007) are used for comparison, but be- 
comes comp atible with the obse rvations if the data from other 
authors (here ICafi'au et al.ll2005h are added. Clearly, more data 
on S abundances in metal-poor stars are needed to reduce the 
current uncertainties in GCE modelling of S evolution, and to 
rule out the set(s) of stellar yields which produce results at vari- 
ance with the observations. 



Calcium The production of calcium is directly attributable 
to inc omplete silicon and oxygen burning d Woosley & Weayep 
Il995h . Departures from non-LTE in the line formation for both 
neutral and singly -ionized calcium have been recently investi- 
gated by Mashon kina et al.l (l2007h . for a wide range of metal- 
licities of the stars. As for Cal, non-LTE abundance corrections 
depend on Teff , log g, [Ca/H] and microturbulence value, and dif- 
fer in value and sign for diff'erent Ca I lines. For Ca II, non-LTE 
leads to negative abundance coiTections over the whole range of 
stellar parameters. 

In Fig. [14] panel c, our model predictions on the behaviour of 
[Ca/Fe] as a function of [Fe/H] are compared to LTE abundances 
from [Cayrel et al. (2004, stars), Lai et al. (2008, pentagons) , 
[bratton et alJ (1200 3. upside -down triang les), lReddv et al i (120031 
[2OO 6, squares) and iBensbv et al.l (l2005 l circles). The uncertain- 
ties in the model predictions for Ca are of the same order of mag- 
nitude as those for Si and S. At very low metallicities. Models 1 
and 5 [short-dashed (red) line and lower solid (black) line, re- 
spectively] produce the best-fitting curves, whereas Model 4 [up- 
per solid (blue) curve] overestimates the observed [Ca/Fe] ratios. 
None of the models is able to properly fit the thin-disc data. 



Titanium The main isotope of titanium, "*^Ti, is produced 
mainly during complete and incomplete silicon burning. 
Spectroscopists may rely on several lines for deriving Ti abun- 
dances in stars. To our knowledge there is only one pub- 
lished non-LTE study of Ti in cool dwarfs and giants by 
iHauschildt et al.l (Il997l) : the data shown in Fig. [14] panel d, are 
all derived under the assumption of LTE. 

As for K (see Sect. I4.2l i. also for Ti the theoretical predic- 
tions are offset with respect to the location of the data points in 
the [Ti/Fe] versus [Fe/H] diagram, though to a lesser extent (see 
Fig. [14] panel d). This clearly points to the need for a major re- 
vision of current SNII and HN yields, over the whole range of 
initial metallicities of the progenitor stars. 



4.4. Iron-group elements 



Scandium Scandium is produced in the innermost ejected lay- 
ers of core-collapse SNe. It is made both as "^^Sc during neon 
burning and as the radioa ctive progenitor "^^Ti in exp losive oxy- 
gen and silicon burning (IWooslev & Weave^ [ 19951) . It IS now 
widely recognized that artificially induced explosions in the pro- 
genitor star models lead to incoiTect nucleosynthesis predictions, 
a consideration which more gene rally apphes t o the composi- 
tion of the whole innermost ejecta dFrohlich et al..,2006> and ref- 
erences therein). A quantity indispensable to correctly describe 
the nucleosynthesis in the innermost ejecta is the electron frac- 
tion in the layers undergoing explosive Si burning, Y^. A value 
of in excess of 0.5 is obtained with a consistent treatment of 
neutrino inte ractions. This b rings about large amounts of Sc, Ti 
and Zn (e.g. iFrohlich et alj |200d). thus solving the problem of 
their underproduction of the values observed in low-metallicity 
stars. Alternatively, an enhanced Sc p roduction is obtained with 
a low-density model dUmeda & Nom oto 2005), where a weak 
jet expands the interior of the core-collapse SN progenitor be- 
fore a strong jet develops, which explodes the star through a 
strong shock. Ti, Mg, Ca, Co and Zn production is enhanced 
in this model as well. The abundances of Sc, as well as Mn 
and Co, in the SN ejecta are then further enhanced by neutrino- 
nucleu s interactions (the v -process) occuiTing in Si-burning re- 
gions (lYoshida et al.ll20()8h . Unfortunately, these studies do not 
provide complete grids of stellar yields useful for GCE models. 

As a matter of fact, the predicted trend of [Sc/Fe] as a func- 
tion of [Fe/H] severely disagrees with the [Sc/Fe] ratios mea- 
sured in solar neighbourhood stars if the Kobavashi et al.' (2006)) 
yields for massive stars are used, independently of the assumed 
HN fraction (solid lines in Fig. [15] panel a). Marginal agree- 
ment between model predictions and observations is obtained by 
adopting the metallicity-dependent yields of Woosle v & Weaver! 
(II995L case B) for massive stars, but only for [Fe/H] < -1.5 
[short-dashed (red) line in Fig. [15] panel a]. 



Vanadium Vanadium is synthesiz ed mainly by incornplete ex- 
plosive Si burning in massive stars dLimongi & Chieffill2003h . 

The difficulty in measuring the V abundances in very metal- 
poor stars makes it troublesome to try and determine the trend 
of [V/Fe] versus [Fe/H] below [Fe/H] a -1.5. It is clear though 
that current stellar yields underestimate the overall V abundance 
in the Galaxy, by at least 0.3-0.4 dex at disc metallicities (see 
Fig. [15] panel b). 



Chromium Chro mium is made chiefly du ring in complete ex- 
plosi ve Si burning dWooslev & Weaver 1995: Limon gi & Chieffil 
|20 03i). Several autho rs (e.g. iHonda et al...2004: Lai et al.ll20()8F 
Bo nifacio et a n i2009l) have noticed that an offset exists between 
the abundances derived from CrI and CrII lines. Furthermore, 
an offset is found also between the abundances measured in gi- 
ants and turnoff stars. The discrepancies may indicate non-LTE 
effects, unaccounted for in current analyses. 

In Fig. \T5\ panel c, we display the predicted behaviour of 
[Cr/Fe] versus [Fe/H] (curves), together with the data from dif- 
ferent authors (symbols). Model 5 [lower solid (black) line], in- 
cluding HN nucleosynthesis, is the one in better agreement with 
the observations. It favours the abundances derived from CrII 
lines, rather th an from Cr I, as indicat ive of the Cr content of halo 
stars (see also IKobavashi et al]|2006h . Notice that CrII lines can 
only be observed in giants. 
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Fig. 15. Predicted behaviour of the relative ratios of several Fe-peak elements to iron (panel a: scandium, panel b: vanadium, 
panel c: chromium, panel d: manganese, panel e: cobalt, panel f: nickel) with respect to the relative iron abundance in the solar 
neighborhood, for Models 1 [short-dashed (red) lines], 4 [upper solid (blue) lines for Sc, V, Cr and Mn; lower solid (blue) lines for 
Co a nd Nil and 5 Flowe r solid (black) lines for Sc, V, Cr and Mn ; upper solid (black) lines for Co and Ni]. Also shown are data 
from Cavrel et al. (2004, stars). iLai et al ] 42008, open pentagons), Bonifacio et al. (2009, asterisks, for Cr from Cril lines, see text), 
Bergemann & Gehren (20081 open circles), Bergemann et al. (2010, open triangles), Gratton et al. (2003, upside-down triangles), 
Reddv et al. (2006, filled squares). iReddv et al.l ( l2003i open squares), i Bensbv et al. (200 5), filled circles, onl y for Cr and Ni) and 
Feltzing et al.. (,2007L filled circles, only for Mnt. i Bergernann & Gehren (^008) and Bergemann et al.l (1201 Oh performed non-LTE 
calculations; all the other abundance studies, instead, rely on the assumption of LTE. See the electronic edition of the journal for a 
colour version of this figure. 



Manganese It is well established through theoretical modelling 
that manganese is mainly made i n explosive s ihcon burning 
and nuclear statistical burning (e.g. ^ Arnettl[r99 6'). but its yields 
from different types of SNe remain still uncertain. It is also 
debated whether or not the yields are m etallicity dependent 
(iFeltzing et al.ll2007h . fCescutti et all (l2008i) have demonstrated 
that the evolution of Mn can be understood only if Mn yields 



from SNela are metallicity-dependent. This is a robust result, be- 
cause it explains at the same time the evolution of Mn in the solar 
vicinity, the Galactic bulge and the Sagittarius dwarf spheroidal 
galaxy. 

The Mn line at 602.1 nm is poorly reproduced by available 
hyper fine structure (HFS) lin elist and should henc e be discarded 
in careful abundance studies (IFeltzing et al.ll20()7h . iReddv et al.l 
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include this line in their analysis, which may explain why 
they find similar trends of [Mn/Fe] versus [Fe/H] for their thick- 
disc and thin-disc samples, while Feltzing et al. (2007) find dif- 
ferent trends below [Fe/H] = (stars with kinematics typical of 
the thick disc show a steady increase of [Mn/Fe] with increasing 
[Fe/H], while stars with kinematics typical of the thin disc share 
almost the same [Mn/Fe] ratio). To our knowledge, there is only 
one work dealin g with the important issue of non-LTE line for- 
mation for Mn dBergemann & Gehrenll2008l) . Though only 14 
stars with [Fe/H] mostly from -1 to -2.5 are analyzed, there 
is intriguing circumstantial evidence that non-LTE corrections 
might run up to 0.5-0.7 dex at low metallicity (Fig.[T5j panel d, 
open circles). 

Let us now compare in Fig. [15] panel d, our model pre- 
dictions on the behaviour of [Mn/Fe] as a function of [Fe/H] 
(solid and dashed lines) to the data obtaine d under the as- 
sumption of LTE by Cavrel et al. (2004, stars). iLai etall (l2008l 
open pentag o ns). iGratton et al.. (,2003^ upside-down triar igles), 
iReddv et al.1 (I2003L l2006l "squares) and iFeltzing et al] (l2007l 
filled circles). At the low metallicities typical of halo stars, say 
[Fe/H] < -1.5, the uncertainties in the GCE model predictions 
are very large, up to 0.8 dex. A model intermediate between 
Models 4 and 5 [upper (blue) and lower (black) solid lines, re- 
spectively], i.e. one allowing a significant fraction of high-mass 
stars to explode as HNe, can explain the observed [Mn/Fe] ra- 
tios very well. However, whe n considering the non-LTE c orrec- 
tions to the Mn abundances dBergemann & Gehrenll2008l open 
circles), the models without HNe (Models 1 and 4) better fit the 
data. At higher metallicities, and especially for [Fe/H] > -0.8, 
all the models fail to reproduce the observations. This is due 
to the abrupt rise of the [Mn/Fe] ratio predicted owing to the 
late, overwhelming contribution to Mn synthesis from SNela. 
Compelling evidence for a reduction of the global Mn yield 
from SNela is discussed in Cescutti et al. (2008, and references 
therein). However, we feel that full non-LTE analysis of larger 
samples of stars is badly needed before we can draw any firm 
conclusion about the chemical evolution of Mn. 

Cobalt Cobalt is synthesized by complete Si burning in the 
deepest stellar layers, whereas Cr and Mn form in the outer in- 
complete Si-burning regions. Up to a few years ago, it was com- 
mon wisdom that the [Co/Fe] ratio as measured in halo stars 
steady increases with decreasing metallicity for -4.0 < [Fe/H] < 
-2.5, whereas the [Cr/Fe] and [Mn/Fe] ratios decrease with de- 
creasing [Fe/H] over the same metallicity range in the same ob- 
jects. The difiference in the observed trends wa s first explained 
by the mass-dependent mass cut in SNell by iNakamura et al.l 
(119991) . Then, Umeda & Nomoto (2005) demonstrated that en- 
hanced explosion energies may lead to higher Co and lower Cr 
and Mn production in low-metallicity, high-mass stars, in agree- 
ment with the observations. However, as mentioned in previous 
paragraphs, the usage of the Cr II lines to estimate Cr abundances 
in giant halo stars and the correction of Mn abundances for non- 
LTE effects have led to a radically diff'erent picture, in which 
Cr and Mn stay almost constant rather than decline towards the 
lowest metallicity end. At the same time, investigations of the 
influence of non-LTE and HFS splitting on the formation of Co 
lines and abundances in cool stars have resulted in a trend of 
increasing [Co /Fe] with decreasing [F e/H] much steeper than 
thought before dBergemann et al.l2010 , open triangles in Fig.fTSl 
panel e). 

We have tried different prescriptions on the synthesis of Co 
in stars in the framework of our GCE model, but none of them 



has brought the model predictions into agreement with the rele- 
vant data (Fig. [15] panel e, solid and dashed lines). This is spe- 
cially true if the non-LTE corrected values are used for compar- 
ison. 

Nickel Nickel is mainly produced in the zo nes which undergo 
complete explosive Si burning in stars (e.g. iLimongi & Chieffj 
2003). Measurements of Ni in solar neighbourhood stars indicate 
a flat trend of [Ni/Fe] versus [Fe/H] over the whole metallicity 
range (Fig. [15] panel f, symbols). 

None of our GCE models is able to explain the observed 
data (Fig. [15] panel f, solid and dashed curves). The uncer- 
tainties in the GCE model predictions are particularly severe at 
the lowest metallicities, where the signature of SNII nucleosyn- 
thesis clearly shows up. Partial agreement between theoretical 
predictions and observations is found with Model 5, adopting 
iKobayashi et al. (2006) yields for massive stars with e = 1, but 
only in the low metallicity domain. The jump in [Ni/Fe] seen 
at [Fe/H ] ~ ~1 is due to Ni production from SNela (see also 
Kobavas hi et al.l l2006l) . The Ni yield from SNela is a strong 
function of 7^ in the burning region, which, in turn, is deter- 
mined by electron captures and, hence, is sensitive to uncertain 
quantities such as the propagation speed of the bu rning front and 
the central density of the white dwarf progenitor dlwamoto et al.l 
1999). This leaves room for a downward revision of current Ni 
yields from SNela dKobavashi et al.ll2006h . 

Copper iRomano & Matteuccil d2007l) have recently presented 
a comprehensive study of copper evolution in different systems. 
They have demonstrated that, on a galactic scale, Cu is mainly 
produce d by massive stars, rather than b y SNela as previously 
thought dMatteucci et al.ll 19931: iMishenina et al.ll2002l) . In mas- 
sive stars, Cu is made during core-helium and carbon-shell hy- 
dros tatic burnings, as well as in explosive comple te Ne burning 
(e.g. lWooslev & Weaveiill99.5HLimongi & Chie fli''2003). 

In Fig. [16] panel a, we compare the predictions from 
Models 1 [short-dashed (red) line], 4 [solid (blue) line labelled 4] 
and 5 [solid (black) line labelled 5] to Cu data from a number of 
studies spanning the whole metallicity range from [Fe/H] ^ -4 
to [Fe/H] ^ -1-0.5. In low-metallicity stars, Cu abundances are 
derived under the assumption of LTE from the near-UV lines of 
Cu I at 3247.53 and 3273.95 A, taking the effect s ofHFSandiso - 
topic splitting into account (Bihain et al. 2004; Lai et al. 2008^. 
At higher metallicities, Cu abundances (Reddy et al. 2003, 200q) 
are obtained from LTE analysis of optical Cul lines (5105.550, 
5218.210 and 5220.090 A); thorough consideration of the HFS 
and isotopic sp litting is required for the 5105 A line. As already 
pointed out by iRomano & Matteuccil d2007h . the Cu content of 
extremely metal-poor stars is basically determined by explosive 
nucleosynthesis in massive stars, but as the metallicity grows, 
the contribution from the weak i-process operating in massive 
stars becomes increasingly important. 

The uncertainties in modelling the Galactic evolution of Cu 
range from 0.9 dex at the lowest metallicities to 0.3 dex or 
e ven less above [Fe/H] = -2.5. Apparently, Model 1, adopting 
iWooslev & Weavei] dl995h case B yields for normal SNell, pro- 
vides the best fit to the observed [Cu/Fe] versus [Fe/H] trend 
for [Fe/H] > -2. At lower metallicities, a model intermediate 
between Models 4 and 5, accounting for a fraction of HNe, re- 
produces better the data. However, the reader should be aware 
that Cu formation through the main i-process acting in low-mass 
AGB stars has not been included in the models, following the 
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results of s tellar nucleosynthesis studie s that found it to be neg- 
ligible (see Romano & Matteucci 2007, and references therein). 
We feel that, before drawing any firm conclusion on which yield 
set would be better for Cu, it would be worth re-investigating 
the issue of Cu production in LIMSs by means of updated AGB 
models. Low-metallicity intermediate-mass AGB stars could in 
principle produce enough Cu to affect the chemical evolution of 
this element. For example, a 5-Mq AGB star with [Fe/H] w -2.3 
produces [Cu/Fe] ~ 1 .0 dex. Cu is produced by neutrons released 
during convective thermal pulses by the ^-Ne(a, n)^^Mg reac- 
tion, similar to the weak ^-process in massive stars. At present, 
yields of Cu from a full range of AGB masses and metallicities 
are not available to test this with a chemical evolution model. 

Zinc Similarly to sulphur, zinc is seen in damped Lyman a sys- 
tems where, due to its refra ctory beh aviour, constitutes a sen- 
sible nucleosynthetic tracer (iNissen e t al. 2004, and references 
therein). Observations of stars in the solar vicinity show that 
[Zn/Fe] is close to ze ro for metallicities in the interval -2.0 < 
[Fe/H] < -1-0.4 (e.g. iBensbv et all l2005t iNissen et all l2007l) . 



but rises steeply to [Zn/Fe] ^ +0.5 at the lowest metallicities 
(ICavreletalJl2004 . Such a trend is reinforced when non-LTE 
corrections are applied (Takeda et al. 2005). An extreme alpha- 
rich freezeout with > 0.5 is a practical mech anism to make 
zinc in high abundance from metal-poor SNell (iFrohhch et al.l 
2006). An alternative possibility for high Zn/Fe production are 
Pop III core-collapse very massive {m ~ 500-1000 M0) stars, 
in which, at the presupernova stage, silicon-burning regions oc- 
cupy a large fractio n, more than 20 per cent of the total mass 
(IOhkuboetal.ll2006h . 

In Fig. [16l panel b, we compare the predictions from 
Models 1 [short-dashed (red) line], 4 and 5 (solid lines labelled 
4 and 5, respectively) to LTE Zn abundances from several stud- 
ies (symbols; see figure caption). Large uncertainties are asso- 
ciated to the model predictions at metallicities typical of the 
halo phase. None of the models is able to reproduce the obser- 
vations: Model 1 is in qualitative agreement with the observed 
trend of increasing [Zn/Fe] with decreasing [Fe/H] in the metal- 
licity interval -4 < [Fe/H] < -1, but it underpredicts Zn. In con- 
trast, the [Zn/Fe] versus [Fe/H] relation predicted by Model 5 for 
[Fe/H] < - L5 is flat. Finally, Model 4 produces by far less Zn 
than required to explain the observations of this element in solar 
neighbourhood stars. 
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Fig. 16. [Cu/Fe] versus [Fe/H] (upper panel) and [Zn/Fe] versus 
[Fe/H] (lower panel) relations as predicted by Models 1 [short- 
dashed (red) lines], 4 and 5 (solid lines labelled 4 and 5, respec- 
tively) for the solar neighb ourhood. Also sh own are data from 
Cavrel et aD (l2004l stars), iLai et al. (l2008l open pentagons). 
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crosses), Bihain et al. (l2004t open trian 



gles), [G ratton et al. ( 2003, upside-do wn triangles), iReddv et aP 
(2006, filled squares), iReddv et all (12003 '. open squares) and 
iBensbv et al.. (2005^. filled circles). See text for discussion. 



5. The Galactic gradient 

The distribution with Galactocentric distance of the abundances 
predicted by a GCE model depends not only on the adopted 
yields, but even more on the assumptions on the galactic param- 
eters, star formation and infall laws and IMF. In particular, the 
slope of the gradients is strictly related to the radial variation 
of the star forrnation/infall rates, i.e. of the enrich ment/dilution 
rates (e.g. iTosil l 19881: iMatteucci & FranqoislllQSQl) . The yields, 
however, provide the absolute abundances and are therefore im- 
portant ingredients for the gradients too. 

The intriguing theme of the establishment and evolution of 
abundance gradients will be the subject of the next paper in this 
series. Here, for the sake of completeness we simply comment 
on the predicted present-day gradient for a few elements of in- 
terest. In Fig. [171 we show the present-day gradients of carbon, 
oxygen, sodium, magnesium, aluminium and silicon predicted 
by Models 1, 3, 4, 5, 9 and 15. For the comparison with the 
observations, we use the homogeneous data for nearly 200 clas- 
sical Cepheid variable stars presented in a series of papers by 
S. Andrieysky and collaborators (lAn drievskv et al. 2002a b,d; 
Luck et al."2003'; ' Andrievskv e"taril2004i: iKovtvukh et al.»2005l: 
Luck et al. 2006). 

It is clear that for some elements different choices of the 
yields do not provide only different absolute abundances, but 
also different slopes of the gradients. This is the case for carbon 
and oxygen, while a negligible effect is seen for the other ele- 
ments, including elements heavier than Si, not shown in Fig. [17] 
It is worth noticing, in particular, that the models including mass 
loss and rotation in high-mass stars are the ones that better re- 
produce the slope of the present-day C gradient. As an example, 
in Fig. [17] panel a, we show the predictions of Model 9, using 
the Geneva group yields for rotating massive stars (upper dot- 
ted line). The theoretical curve is arbitrarily shifted by -0.5 dex 
(lower dotted line) to better compare the slope of the theoretical 
gradient with that of the gradient observed in Cepheids. Other 
models, using different nucleosynthesis prescriptions, predict a 
much flatter C gradient, at vari ance with observations. Model 3, 
using the older Maeder ( Il992h yields for solar-metallicity mas- 
sive stars computed with high mass loss but no rotation, is in 
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Rq [kpc] R3 [kpc] 

Fig. 17. Radial distribution of the abundances of carbon, oxygen, sodium, magnesium, aluminium and silicon in the disc of the 
Milky Way at the present time. The predictions from Models 1 [short-dashed (red) lines], 3 [long-dashed (black) lines, only for C 
and O], 4 [upper solid (blue) curves], 5 [lower solid (black) curves], 9 [dotted (black) lines, only f or C, O and Na] and 15 (t hick 
solid lines, only f o r Na and Al) are compar e d to the homogeneous d ata f or Cepheid variable stars bv lAndrievskv et al] (I2002allbll3) . 
iLuck et al.l (l2003h . lAndrievskv et a"D (l2004 . Ikovtvukh et dl ('2005') and 'Lu ck et"an (|2006V In panels a and d the predictions from 
Model 5 are not shown because they are pretty much the same as Model 4. The intersection of the short-dashed vertical and 
horizontal (grey) lines in each panel marks the position of the Sun. See text for details. See the electronic edition of the journal for 
a colour version of this figure. 



agreement with the slope of the observed C gradient, but pro- 
duces a too flat O gradient [Fig. [17] panels a and b, long-dashed 
(black) lines]. 

There are basically two conditions that the models must meet 
to be judged in agreement with the observations. First, they have 
to pass slightly above the intercept of the dashed lines indicat- 
ing the position of the Sun in each diagram. Second, they have 
to reproduce the slope of the observed gradient. The intersec- 
tion of the vertical and horizontal dashed lines marks the po- 



sition of the Sun in each panel. If the composition of the Sun 
is representative of the typical chemical enrichment occuiTing 
in the Galactic disc at a distance of Rc = ^g.O - 8 kpc from 
the Galactic centre, one expects the abundances of Cepheids at 
comparable Galactocentric distances to scatter slightly above the 
solar value, under the assumption that the metallicity of the ISM 
has only slightly evolved from the time of the Sun formation up 
to now (.Tosi,1988.) . This is what is found for most elements, 
with the notable exceptions of carbon (on average lower than so- 
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Fig. 18. Y versus 10'* O/H in the solar vicinity (left panels) and in the inner disc (Re = 6 kpc; right panels), as predicted by models 
with different nucleosynthesis prescriptions. Top panels: Model 1 [short-dashed (red) lines] versus Model 2 [solid (black) lines], 
showing the effect of changing the mass cut in models of SNll explosions. Bottom panels: Model 6 [dot-dashed (black) lines] 
versus Model 7 [solid (black) lines], with and without rotation in high-mass stars, respectively. Left panels: the filled circles on each 
theoretical track indicate the Y and lO'* O/H values predicted at the time of the birth of the Sun, 4.5 Gyr ago. These are compared 
wit h the bulk (i nitial) abundances of He and O in the Sun (Sun symbols, from Grevesse & Sauval 1998 - higher O abundance - 
and lAsplund et al][2009l - lower O abundance) . The effects of diffus ion are taken into account in defining the protosolar abundances 
displayed in the figure (see the discussion in lAsplund et al.ll2009h . Right panels: the final Y and 10'* O/H values attained by the 
models are highlighted (filled circles) and compared to the abundances of He and O measured for M 17, an H II region located at 
2.1 kpc distance from the Sun (filled triangles; see text for discussion). 



lar in Cepheids, becau se of the occurrence of the first dredge-up; 
iKovtvukh et al.ll2"005h . sodium and aluminium (on average en- 
hanced with respect to solar, again reflecting evolutionary pro- 
cesses inside the stars themselves; Kovtvukh et al. 2005). While 
this may complicate the comparison with the absolute abun- 
dances of C, Na and Al predicted by the models, our conclusions 
on the slopes of the gradients should remain basically unaffected. 

In conclusion, we have seen that, at least for some elements, 
the comparison of the gradients predicted by using different 
grids of yields with a homogeneous data set can prove useful 
to discriminate among different nucleosynthesis prescriptions. 

6. The helium-to-heavy element enrichment ratio 

In Figs.[T8]to|20]we show the predicted behaviour of the helium 
mass fraction, Y, versus 10"* O/H in the ISM for two Galactic 
radii, Rq = ^g,0 - 8 kpc (left panels) and Rc - 6 kpc (right 
panels). Shown are the predictions of models computed with 
different nucleosynthesis prescriptions (see Figs. [18] to |20] cap- 
tions and detailed description of the adopted nucleosynthesis in 
TableEJ. 



The predictions of GCE models on the composition of the 
ISM in the solar vicinity 4.5 Gyr ago are usually compared, 
for elements heavier than He, to the photospheric solar abun- 
dances. However, this is incorrect, since the photospheric abun- 
dances are no longer representative of the original composi- 
tion of the Sun at birth , because of the effects of diffusion 
(e.g. iBahcall et al.l l2006h . The solar abundances of He and O 
displayed in Figs. [18] to [20] for comparison with the theoret- 
ical model results (Sun symbols) have been corrected for the 
effects of diffusion ( Asplund et al. 2009). W e show both the 
values recommended by iGrevesse & Sauval (Il998h . Yq = 
0.2701, log£(O)0,i„i = 8.87 [log g(0)fr, phn, = 8.83 + 0.06], and 
those recommended by lAsplund et"aLr(l2009h . Yq = 0.2703, 
log£(0)o,in, = 8.73 [log e(0)o, phot = 8.69 ± 0.05]' 

It can be immediately seen that the largest uncertainties in 
the model predictions are associated to the predicted oxygen 
abundances, while all the models predict virtually the same Y 
value at /?G = 8 kpc at Sun's birth and at Rc -6 kpc at the 
present time. Small differences are found in the model predic- 
tions when changing the prescriptions on He and O synthesis 
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Fig. 19. Same as Fig. [18] but confronting the predictions from Models 6, 8 and 9 (upper panels; dot-dashed, solid and dotted lines, 
respectively) and 10 and 1 1 (lower panels; short-dashed-long-dashed and solid lines, respectively). The different model predictions 
overlap. This suggests that the uncertainties associated with the treatment of mass loss and HBB in AGB do not have an impact on 
the predicted Y versus lO'^O/H behaviour 



in LIMSs (Fig. [19] and lower panels of Fig. l20b . On the other 
hand, the differences become significant when changing the pre- 
scriptions on He and O synthesis in massive stars (Fig. [TS] and 
upper pa nels of Fig. [20|. In p articular, the models adopting the 
yields bv lHirschi et aL " oOOsI) for solar-metallicity rotating mas- 
sive stars (Models 6, 8, 9, 10, 11, 13 and 15) tend to overpro- 
duce oxygen in the Galaxy. In particular, the helium-to-metal 
enrichment ratio predicted by Model 1 3 around and above so- 
lar metallicity is AF/AZ ^ 1. This is barely consiste nt with 
the va lue of AF/AZ = 2.1 + 0.9 suggested by Casagran de et alj 
(I2007h for a large sample of nearby K dwarf stars. Models 3, 12 
and 14 (t he latter not sh own in the figures), adopting the older 
yields bv iMaederl (Il992h for solar-metallicity stars, instead, are 
consistent with Yf-, = 0.270 3, loge(O)0 - 8.73, as sug- 
gested by Asplu nd et alJ (l2009l) . and give AF/AZ = 1.3-1.5 for 
solar- and supers ol ar-met allicity stars, in good agreement with 
ICasagrande et alJ s (l2007h findings. 



Th e best way to discriminate between the i Hirschi et alJ 
(l2005h and Maeder (1992) yields for solar-metallicity stars 
is to compute the chemical evolution at inner radii where, 
owing to the higher metallicities attained, the differences in 
the model predictions are expected to exacerbate. In the in- 
ner Galaxy, at a distance of 2.1 ± 0.2 kpc from the Sun 
(iHofFmeister e t al. 2008), lies M 17, the best Galactic H 11 region 
to derive He and O abundances, owing to t he smallest correc- 
tion for the presence of neutral helium (see ICarigi & Peimberi 



I2OO8I, and references therein). In Fig. [20] upper right panel, 
we compare the predictions from Models 1 and 3 for Rq = 
6 kpc with t he new determina t ion of He and O abundances 
for M 17 b y ICarigi & PeimbertI (|2008|) . It is seen that the use 
of iMaedeil i[l992h yields for solar-metallicity massive stars im- 
proves the model predictions. Independent chemical evolution 
studies have already supported Maeder ( 19 ^ yie l ds for solar- 
metallicity massive star s (Giovagnoli & Tqs| 1995t Carigi et aP 
120051: ICarigi & Peimberti,2008.: .McWilham et al.l l2008f). On the 
other hand, it is worth noticing that also with the up-to-date 
yields computed by the Geneva group for non-rotating massive 
stars, the fits to the Y, O/H values at Rq - ^ kpc at Sun's birth 
and at /?G = 6 kpc at the present time improve with respect to the 
case with rotation (see Fig. [18] lower panels). 

We deem it necessary to stress at this point that different 
assumptions about the model ingredients can change the con- 
clusions. As an example, in Fig. [21] we show the results of 
Models IS and lOS computed for Rc - 6 kpc (top panel) and 
Rg - 8 kpc (bottom panel), this time assuming a Scalo (198^ 
IMF. Because of the highest mass fraction per stellar genera- 
tion locked up in stars that do not contribute to the chemical 
enrichment and because of the lower mass fraction falling in the 
8-100 M0 mass range according to this IMF, the models now 
produce less metals (and, hence, less oxygen) and less helium. 
Thus, while still consistent, within the uncertainties, with the ini- 
tial He and O abundances in the Sun, they are by no means able 
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Fig. 20. Same as Fig.[T8j but confronting the predictions from Models 1 and 3 [upper panels; short-dashed (red) and long-dashed 
(black) lines, respectively] and 13 and 15 (lower panels; thick dotted and thick solid lines, respectively). 



to account for the abundances of He and O in M 17. Adopting 
different nucleosynthesis prescriptions would only worsen the 
discrepancy between model predictions and observations. 

Finally, it is worth noticing that no constraints can be 
set at present on the helium-to-metal enrichment ratio in 
the solar neig hbourhood at ages older than that of the Sun 
dCasagrande e t al. 2007). This means that no constraints can be 
set on the He stellar yields for Z < Zq from chemical evolution 
studies of the Galaxy. 



7. Summary and conclusions 

In this paper, we have described the evolution of He and heavier 
elements from C to Zn in the Milky Way resulting from different 
sets of stellar yields and a well-tested chemical evolution model 
for the Galaxy. 

In Sect. |2] we have provided a detailed comparison among 
different yield sets suitable for use in chemical evolution models 
of the Galaxy. We stress that the absence of a complete and ho- 
mogeneous grid of stellar yields, computed with the same input 
physics for many initial masses and chemical compositions of 
the stars - say, from 0.8 Mq to 100 M0 and from zero to super- 
solar metallicity - is a major cause of uncertainty in the model 
predictions for elements such as He and N. In particular, detailed 
grids of yields for super- AGB stars are badly needed. 

In Sect. [3] we introduced the adopted model for the chem- 
ical evolution of the Milky Way. This model well reproduces 
all the main observed features of the solar vicinity and of the 



whole Galactic disc. In particular, it is in agreement with the cur- 
rent star formation rate, gas and star radial density profiles and 
with the present-day Galactic SNll and SNla rates and their ra- 
tio; it also reproduces the relative number of metal-poor-to-total 
stars, infall rate and degree of deuterium astration in the solar 
neighbourhood (see IChiap pini et al.l 120011; iRomano et al.l 120061 
and Paper 1). All these quantities are reproduced independently 
of the choice of the stellar yields and with a minimum number of 
free parameters, which justifies the use of the code to test possi- 
ble stellar nucleosynthesis scenarios. 

In Sects. |4]to|6] we described the model results, obtained by 
assuming different sets of stellar yields. In particular, in Sect. |4] 
we dealt with the behaviour of several abundance ratios as func- 
tions of metallicity in the solar neighbourhood. The radial abun- 
dance profiles were discussed in Sect. |5] whereas in Sect. |6] 
we focused on the relative helium-to-metal enrichment in the 
Galaxy. A meaningful comparison between model predictions 
and observations relies on carefully selected, reliable data sets. 
Hence, in Sects. |4]to|6] we also tried to ascertain the origin of 
possible inconsistencies on different data sets. We find that, in 
many cases, non-LTE and/or 3D analyses of stellar spectra lead 
to a profound revision of the abundance trends inferred by means 
of ID, plane-parallel model atmospheres under the LTE approx- 
imation. This deeply impacts on the interpretation of the model 
results. 

Despite several considerable improvements in the field of 
stellar evolution and nucleosynthesis in recent years, no single 
combination of stellar yields is found which is able to repro- 
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duce at once all the available measurements of chemical abun- 
dances and abundance ratios in the Milky Way. In Figs. l22l 
and l23l we summarize the results of our 'best yield choice', 
i.e. the combination of yields which maximizes the agreement 
with the largest number of observational constraints for the 
Milky Way. Our best model. Model 15 (solid lines in Figs. 1221- 
1231 ). uses the yields by Karakas (2010) for LIMSs, the pre- 
supernova yields of He and CNO elements by the Geneva group 
for massive stars , and the results of explosive nucleosynthesis of 
iKobavashi etall (|2006) for SNell below 20 Mq and HNe above 
that limit. For sake of comparison, we also show the results 
of Mode l 1 (dash ed lines i n Figs. [22l423l l. computed with the 
yields of van den H oek & Groenewegen ( 1997) for LIMSs and 
IWoosIev & Weavei] (119951) for massive stars. This is the combi- 
nation of stellar yields most often used in the literature. 

Model 15 fits better the behaviour of [C/Fe], [N/Fe], [Na/Fe], 
[Mg/Fe], [Al/Fe], [Cr/Fe] and, perhaps, [Zn/Fe] as a function of 
[Fe/H] in the solar neighbourhood, although some caveats still 
remain: 

1. An exceptionally good fit between model predictions and 
observations of carbon in the disc is obtained if the car- 



bon abundances are those based on the forbidden [C I] line 
at 872.7 nm (Bensbv & Feltzing 2006, golden points in 
Fig.l22l): otherwise, the fit is rather poor 

2. The overall trend of [N/Fe] versus [Fe/H] can not be satis- 
factorily explained by any model yet. 

3. The behaviours of Na in halo stars and of Mg and Al at disc 
metallicities are still far from being understood. Rotation 
might strongly impact on yields of Na and Al from mas- 
sive stars, but detailed computations suitable for use in GCE 
models of the Galaxy are still missing in the literature. 

The evolution of oxygen in the solar neighbourhood and its so- 
lar abundance are very well reproduced by both models, thus 
indicating that its nucleosynthesis in stars is well understood by 
now. The same is true for Si, while S and Ca are slightly under- 
estimated in the disc. 

Potassium, scandium, titanium and vanadium are understi- 
mated by both models, with Model 15 performing worse for K, 
Sc and Ti. Ni is overestimated by Model 1, whereas Model 15 
partially agrees with the data for halo stars. 

Should future non-LTE analy ses of larger samples of stars 
confirm the trends suggested by iBergemann & GehrenI (l2008h 
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Fig. 23. [X/Fe] predicted by Model 1 (filled circles and dashed line) and Model 15 (filled squares and solid line) at Sun's birth 
for elements from C to Zn. The horizontal dotted line indicates the solar value of the ratios that the models must reproduce. The 
horizontal dashed Unes allow for a ±0. 15 dex error in the model predictions. 



and iBergemann et alj ( 1201 Ol) for Mn and Co, respectively, we 
would be led to conclude that: (i) none of the models is able to 
explain the Co data in the Galaxy; (ii) Model 1 reproduces the 
halo data for Mn better than Model 15; (iii) neither model can 
explain the behaviour of Mn in the Galactic disc, most likely be- 
cause of an excessive synthe sis of Mn from SNel a at late times. 
Indeed, as demonstrated by ICescutti et alj (l2008l) . who studied 
Mn evolution in the solar neighbourhood, the Galactic bulge and 
the Sagittarius dwarf spheroidal galaxy, considering metallicity- 
dependent yields of Mn from SNela could provide a solution to 
the problem of the overproduction of Mn at late times. 

Finally, Model 1 provides a very good fit to copper abun- 
dances measured in solar neighbourhood stars, while Model 15 
fails to reproduce Cu abundances in disc stars. However, the role 
of AGB stars as Cu factories has still to be assessed, which may 
lead to improvements in Model 15 predictions for Cu at disc 
metallicities. 

While exploding core-collapse SNe through an artificially 
induced piston or thermal bomb is acceptable for the external 
layers, only a coiTect treatment of the physics of the explosion 
could provide coiTect yields for the species produced in the in- 
nermost ejected layers, such as the Fe-peak nuclei. Our results 
are in agreement with recent claims that a correct treatment of 
the explosion physics is still missing in SN models. 

The inclusion of rotation in high-mass stellar models has al- 
lowed to sensibly improve the predictions of GCE models on 
C evolution in the early Galaxy. Contrary to common wisdom, 
high N/O ratios at the lowest metallicities, instead, can be ob- 
tained even if not allowing for low-metallicity massive stars to 
rotate fast. 

Yields for LIMSs computed by taking rotation into account 
are provided in the literature, but stopping the computations at 
the beginning of the TP- AGB stage makes them useless in GCE 



models, since the important effects of the TDU and HBB on the 
final yields are not considered. 

We sincerely hope that the results of this study can trigger a 
further effort to improve and expand the current existing grids of 
stellar yields. 
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